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Abstract. These are introductory lectures to some aspects of the physics of strongly correlated 
electron systems. I first explain the main reasons for strong correlations in several classes of 
materials. The basic principles of dynamical mean-field theory (DMFT) are then briefly reviewed. 
I emphasize the formal analogies with classical mean-field theory and density functional theory, 
through the construction of free-energy functionals of a local observable. I review the application 
of DMFT to the Mott transition, and compare to recent spectroscopy and transport experiments. 
The key role of the quasiparticle coherence scale, and of transfers of spectral weight between 
low- and intermediate or high energies is emphasized. Above this scale, correlated metals enter 
an incoherent regime with unusual transport properties. The recent combinations of DMFT with 
electronic structure methods are also discussed, and illustrated by some applications to transition 
metal oxides and f-electron materials. 
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1. INTRODUCTION: WHY STRONG CORRELATIONS ? 
1.1. Hesitant electrons: delocalised waves or localised particles ? 

The physical properties of electrons in many solids can be described, to a good ap- 
proximation, by assuming an independent particle picture. This is particularly successful 
when one deals with broad energy bands, associated with a large value of the kinetic en- 
ergy. In such cases, the (valence) electrons are highly itinerant: they are delocalised over 
the entire solid. The typical time spent near a specific atom in the crystal lattice is very 
short. In such a situation, valence electrons are well described using a wave-like picture, 
in which individual wavefunctions are calculated from an effective one-electron periodic 
potential. 

For some materials however, this physical picture suffers from severe limitations and 

may fail altogether. This happens when valence electrons spend a larger time around 
a given atom in the crystal lattice, and hence have a tendency towards localisation. In 
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such cases, electrons tend to "see each other" and the effects of statistical correlations 
between the motions of individual electrons become important. An independent particle 
description will not be appropriate, particularly at short or intermediate time scales (high 
to intermediate energies). A particle-like picture may in fact be more appropriate than a 
wave-like one over those time scales, involving wavefunctions localised around specific 
atomic sites. Materials in which electronic correlations are significant are generally 
associated with moderate values of the bandwidth (narrow bands). The small kinetic 
energy implies a longer time spent on a given atomic site. It also implies that the ratio 
of the Coulomb repulsion energy between electrons and the available kinetic energy 
becomes larger. As a result delocalising the valence electrons over the whole solid 
may become less favorable energetically. In some extreme cases, the balance may even 
become unfavorable, so that the corresponding electrons will remain localised. In a naive 
picture, these electrons sit on the atoms to which they belong and refuse to move. If this 
happens to all the electrons close to the Fermi level, the solid becomes an insulator. 
This insulator is difficult to understand in the wave-like language: it is not caused by the 
absence of available one-electron states caused by destructive interference in k-space, 
resulting in a band-gap, as in conventional band insulators. It is however very easy to 
understand in real space (thinking of the solid as made of individual atoms pulled closer 
to one another in order to form the crystal lattice). This mechanism was understood long 
ago [1, 2] by Mott (and Peierls), and such insulators are therefore called Mott insulators 
(Sec. 4. In other cases, such as f-electron materials, this electron localisation affects only 
part of the electrons in the solid (e.g the ones corresponding to the f-shell), so that the 
solid remains a (strongly correlated) metal. 

The most interesting situation, which is also the one which is hardest to handle 
theoretically, is when the localised character on short time-scales and the itinerant 
character on long time-scales coexist. In such cases, the electrons "hesitate" between 
being itinerant and being localised. This gives rise to a number of physical phenomena, 
and also results in several possible instabilities of the electron gas which often compete, 
with very small energy differences between them. In order to handle such situations 
theoretically, it is necessary to think both in k-space and in real space, to handle both 
the particle-like and the wave-like character of the electrons and, importantly, to be 
able to describe physical phenomena on intermediate energy scales. For example, one 
needs to explain how long-lived (wave-like) quasiparticles may eventually emerge at low 
energy/temperature in a strongly correlated metal while at higher energy/temperature, 
only incoherent (particle-like) excitations are visible. It is the opinion of the author 
that, in many cases, understanding these intermediate energy scales and the associated 
coherent/incoherent crossover is the key to the intriguing physics often observed in 
correlated metals. In these lectures, we discuss a technique, the dynamical mean-field 
theory (DMFT), which is able to (at least partially) handle this problem. This technique 
has led to significant progress in our understanding of strong correlation physics, and 
allows for a quantitative description of many correlated materials [3, 4]. Extensions and 
generalisations of this technique are currently being developed in order to handle the 
most difficult/mysterious situations which cannot be tamed by the simplest version of 
DMFT. 
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1.2. Bare energy scales 



Localised orbitals and narrow bands. In practice, strongly correlated materials are 
generally associated with partially filled d- or f- shells. Hence, the suspects are materials 
involving: 

• Transition metal elements (particularly from the 3d-shell from Ti to Cu, and to a 

lesser extent 4d from Zr to Ag). 

• Rare earth (4f from Ce to Yb) or actinide elements (5f from Th to Lw) 

To this list, one should also add molecular (organic) conductors with large unit cell 
volumes in which the overlap between molecular orbitals is weak. 

What is so special about d- and f- orbitals (particularly 3d and 4f) ? Consider the 
atomic wavefunctions of the 3d shell in a 3d transition metal atom (e.g Cu). There 
are no atomic wavefunctions with the same value Z = 2 of the angular momentum 
quantum number, but lower principal quantum number n than n — 3 (since one must 
have I < n — \). Hence, the 3d wavefunctions are orthogonal to all the n = 1 and n = l 
orbitals just because of their angular dependence, and the radial part needs not have 
nodes or extend far away from the nucleus. As a result, the 3d-orbital wave functions 
are confined more closely to the nucleus than for s or p states of comparable energy. 
The same argument applies to the 4f shell in rare earths. It also implies that the 4d 
wavefunctions in the 4d transition metals or the 5f ones in actinides will be more 
extended (and hence that these materials are expected to display, on the whole, weaker 
correlation effects than 3d transition metals, or the rare earth, respectively). 

Oversimplified as it may be, these qualitative arguments at least tell us that a key 
energy scale in the problem is the degree of overlap between orbitals on neighbouring 
atomic sites. This will control the bandwidth and the order of magnitude of the kinetic 
energy. A simple estimate of this overlap is the matrix element: 



In the solid, the wavefunction Xi{r — R) should be thought of as a Wannier-like wave 
function centered on atomic site R. In narrow band systems, typical values of the 
bandwith are a few electron-volts. 

Coulomb repulsion and the Hubbard U. Another key parameter is the typical 
strength of the Coulomb repulsion between electrons sitting in the most localized or- 
bitals. The biggest repulsion is associated with electrons with opposite spins occupying 
the same orbital: this is the Hubbard repulsion which we can estimate as: 



In this expression. Us is the interaction between electrons including screening effects by 
other electrons in the solid. Screening is a very large effect: if we were to estimate (2) 
with the unscreened Coulomb interaction f/(r — r') = e^/|r — r'|, we would typically 
obtain values in the range of tens of electron-volts. Instead, the screened value of U 




(1) 



U ~ 



J drdr'\xL{r - R) P f/.(r - r') \xL{r' - R)\^ 



(2) 
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in correlated materials is typically a few electron-volts. This can be comparable to 
the kinetic energy for narrow bandwiths, hence the competition between localised and 
itinerant aspects. Naturally, other matrix elements (e.g between different orbitals, or 
between different sites) are important for a realistic description of materials (see the last 
section of these lectures). 

In fact, a precise description of screening in solids is a rather difficult problem. An 
important point is, again, that this issue crucially depends on energy scale. At very 
low energy, one should observe the fully screened value, of order a few eV 's, while 
at high energies (say, above the plasmon energy in a metal) one should observe the 
unscreened value, tens of eV's. Indeed, the screened effective interaction W{r,r';(o) 
as estimated e.g from the RPA approximation, is a strong function of frequency (see 
e.g Ref. [5, 6] for an ab-initio GW treatment in the case of Nickel). As a result, using 
an energy-independent parametrization of the on-site matrix elements of the Coulomb 
interaction such as (2) can only be appropriate for a description restricted to low- enough 
energies [6]. The Hubbard interaction can only be given a precise meaning in a solid, 
over a large enery range, if it is made energy-dependent. I shall come back to this issue 
in the very last section of these lectures (Sec. 5.5). 

The simplest model hamiltonian. From this discussion, it should be clear that the 
simplest model in which strong correlation physics can be discussed is that of a lattice 
of single-level "atoms", or equivalently of a single tight-binding band (associated with 
Wannier orbitals centered on the sites of the crystal lattice), retaining only the on-site 
interaction term between electrons with opposite spins: 

^ = - XI ^RR'^Rcj^R'cJ + eoJ^«Ra + l^«RT«Ri (3) 
RR',<T R<7 R 

The kinetic energy term is diagonalized in a single-particle basis of Bloch's wavefunc- 
tions: 

^0 = E Ckc+^ckc ; ek = E ^RR' e^<^-^) (4) 

kcT R' 
with e.g for nearest-neighbour hopping on the simple cubic lattice in d-dimensions: 

d 

Ck = —It E cos(^jufl) (5) 

M=l 

In the absence of hopping, we have, at each site, a single atomic level and hence 
four possible quantum states: 1 0) , | t) , | \) and | ti) with energies , eo and U + Ieq, 
respectively. 

Eq. (3) is the famous Hubbard model [7, 8, 9]. It plays in this field the same role than 
that played by the Ising model in statistical mechanics: a laboratory for testing physical 
ideas, and theoretical methods alike. Simplified as it may be, and despite the fact that 
it already has a 40-year old history, we are far from having explored all the physical 
phenomena contained in this model, let alone of being able to reliably calculate with it 
in all parameter ranges ! 
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1.3. Examples of strongly correlated materials 

In this section, I give a few examples of strongly correlated materials. The discus- 
sion emphasizes a few key points but is otherwise very brief. There are many useful 
references related to this section, e.g [10, 11, 2, 12, 13]. 



1.3.1. Transition metals 

In 3d transition metals, the 4s orbitals have lower energy than the 3d and are therefore 
filled first. The 4s orbitals extend much further from the nucleus, and thus overlap 
strongly. This holds the atoms sufficiently far apart so that the d- orbitals have a small 
direct overlap. Nevertheless, d-orbitals extend much further from the nucleus than the 
"core" electrons (corresponding to shells which are deep in energy below the Fermi 
level). As a result, throughout the 3d series of transition metals (and even more so in the 
4d series), d-electrons do have an itinerant character, giving rise to quasiparticle bands. 
That this is the case is already clear from a very basic property of the material, namely 
how the equilibrium unit- cell volume depends on the element as one moves along the 
3d series (Fig. 1). The unit-cell volume has a very characteristic, roughly parabolic, 
dependence. A simple model of a narrow band being gradually filled, introduced long 
ago by Friedel [14] accounts for this parabolic dependence (see also [15, 16]). Because 
the states at the bottom of the band are bonding-like while the states at the top of the 
band are anti-bonding like, the binding energy is maximal (and hence the equilibrium 
volume is minimal) for a half -filled shell. Instead, if the d-electrons were localised we 
would expect little contribution of the d-shell to the cohesive energy of the solid, and the 
equilibrium volume should not vary much along the series. 



Screening is relatively efficient in transition metals because the 3d band is not too far in 

energy from the 4s band. The latter plays the dominant role in screening the Coulomb 
interaction (crudely speaking, one has to consider the following charge transfer process 
between two neighbouring atoms: ZcPAs + ?>(PAs — > ?id^~^As^ + , see e.g [17] 
for further discussion). For all these reasons (the band not being extremely narrow, 
screening being efficient), electron correlations do have important physical effects for 
3d transition metals, but not extreme ones like localisation. Magnetism of these metals 
below the Curie temperature, but also the existence of fluctuating local moments in the 
paramagnetic phase are exemples of such correlation effects. Band structure calculations 
based on DFT-LDA methods overestimate the width of the occupied d-band (by about 
30% in the case of nickel). Some features observed in spectroscopy experiments (such 
as the (in)famous 6 eV satellite in nickel) are also signatures of correlation effects, and 
are not reproduced by standard electronic structure calculations. 
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Lanthanides: Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb Lu 
Actinides: Th Pa U Np Pu Am Cm Bk Cf 
Transition metals: La Hf Ta W Re Os Ir Pt Au 

FIGURE 1. Experimental Wigner-Seitz Radius of Actinides, Lanthanides, and 5d Transition Metals. 
The equilibrium volume of the primitive unit cell is given by V = 4-KR^g/3. Elements that lie on top of 
each other have the same number of valence electrons. The volume of the transition metals has a roughly 
parabolic shape, indicating delocalised 5d electrons. The volumes of the lanthanides remain roughly 
constant, indicating localised 4f electrons. The volumes of the light actinides decrease with increasing 
atomic number, whereas the volumes of the late actinides behaves similarly to that of the lanthanides. 
FromRef. [16] 



1.3.2. Transition metal oxides 

In transition metal compounds (e.g oxides or chalcogenides), the direct overlap be- 
tween d-orbitals is generally so small that d-electrons can only move through hybridisa- 
tion with the ligand atoms (e.g oxygen 2p-bands). For example, in the cubic perovskite 
structure shown on Fig. 2, each transition-metal atom is "encaged" at the center of an 
octahedron made of six oxygen atoms. Hybridisation leads to the formation of bond- 
ing and antibonding orbitals. An important energy scale is the charge-transfer energy 
A = Crf — £p, i.e the energy difference between the average position of the oxygen and 
transition metal bands. When A is large as compared to the overlap integral tp^, the bond- 
ing orbitals have mainly oxygen character and the antibonding ones mainly transition- 
metal character. In this case, the effective metal- to- metal hopping can be estimated as 
teff ~ fpd/^^ ^^'^ therefore quite small. 

The efficiency of screening in transition-metal oxides depends crucially on the relative 
position of the 4s and 3d band. For 3d transition metal monoxides MO with M to the 
right of Vanadium, the 4s level is much higher in energy than 3d, thus leading to poor 
screening and large values of U. This, in addition to the small bandwidth and relatively 
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FIGURE 2. The cubic perovskite structure, e.g of the compound SrVOs . Transition-metal atoms (V) - 
small grey spheres- are at the center of oxygen octahedra (dark spheres). Sr atoms are the larger spheres 
in-between planes. FromRef. [18]. 



large A, leads to dramatic correlation effects, turning the system into a Mott insulator 
(or rather, a charge- transfer insulator, see below), in spite of the incomplete filling of 
the d-band. The Mott phenomenon plays a key role in the physics of transition- metal 
oxides, as discussed in detail later in these lectures (see Sec. 4 and Fig. 6). 

Crystal field splitting. The 5 -fold (10-fold with spin) degeneracy of the d- orbitals 
in the atom is lifted in the solid, due to the influence of the electric field created 
by neighbouring atoms, i.e the ligand oxygen atoms in transition metal oxides. For a 
transition metal ion in an octahedral environment (as in Fig. 2), this results in a three-fold 
group of states {t2g) which is lower in energy and a doublet {eg) higher in energy. Indeed, 
the dxy,dy^,d^ orbitals forming the t2g multiplet do not point towards the ligand atoms, 
in contrast to the states in the doublet {d^2_y2, dT^^i_^i). The latter therefore lead to a 
higher cost in Coulomb repulsion energy. For a crystal with perfect cubic symmetry, the 
t2g and eg multiplets remain exactly degenerate, while a lower symmetry of the crystal 
lattice lifts the degeneracy further. For a tetrahedral environment of the transition-metal 
ion, the opposite situation is found, with t2g higher in energy than gg. In transition metals, 
the energy scale associated with crystal- field splitting is typically much smaller than 
the bandwith. This is not so in transition-metal oxides, for which these considerations 
become essential. In some materials, such as e.g SrVOs and the other d^ oxides studied 
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in Sec. 5.3 of these lectures, the energy bands emerging from the ?2g and eg orbitals form 
two groups of bands well separated in energy. 

Mott- and charge -transfer insulators. There are two important considerations, 
which are responsible for the different physical properties of the "early" (i.e involving 
Ti, V, Cr, ...) and "late" (Ni, Cu) transition- metal oxides: 

- whether the Fermi level falls within the ?2g or eg multiplets, 

- what is the relative position of oxygen (e^) and transition-metal (Ep) levels ? 
For those compounds which correspond to an octahedral environment: 

• In early transition-metal oxides, t2g is partially filled, eg is empty. Hence, the 
hybridisation with ligand is very weak (because ?2g orbitals point away from the 
2p oxygen orbitals). Also, the d-orbitals are much higher in energy than the 2p 
orbitals of oxygen. As a result, the charge-transfer energy A = — gp is large, and 
the bandwidth is small. The local d-d Coulomb repulsion Udd is a smaller scale than 
A but it can be larger than the bandwidth (~ t^^/A): this leads to Mott insulators. 

• For late transition-metal oxides, t2g is completely filled, and the Fermi level lies 
within eg. As a result, hybridisation with the ligand is stronger. Also, because of 
the greater electric charge on the nuclei, the attractive potential is stronger and 
as a result, the Fermi level moves closer to the energy of the 2p ligand orbitals. 
Hence, A is a smaller scale than Udd and controls the energy cost of adding an extra 
electron. When this cost becomes larger than the bandwidth, insulating materials 
are obtained, often called "charge transfer insulators" [19, 20]. The mechanism is 
not qualitatively different than the Mott mechanism, but the insulating gap is set 
by the scale A rather than Udd and separates the oxygen band from a d-band rather 
than a lower and upper Hubbard bands having both d-character. 

The p-d model. The single-band Hubbard model is easily extended in order to take 
into account both transition-metal and oxygen orbitals in a simple modelisation of 
transition-metal oxides. The key terms to be retained are^: 

Hpd = - £ tpd{d^^PR'a + h.c) + edY,f^Ra + ^pT,f^i'a + ^ddY.'^i^'^il (6) 
RR',<7 RcT R'cT R 

to which one may want to add other terms, such as: Coulomb repulsions Upp and Upd or 
direct oxygen-oxygen hoppings tpp. 



1.3.3. f -electrons: rare earths, actinides and their compounds 

A distinctive character of the physics of rare-earth metals (lanthanides) is that the 4f 
electrons tend to be localised rather than itinerant (at ambiant pressure). As a result, the 



For simplicity, the hamiltonian is written in the case where only one d-band is relevant, as e.g for 
cuprates. 
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f-electrons contribute contribute little to the cohesive energy of the solid, and the unit- 
cell volume depends very weakly on the filling of the 4f shell (Fig. 1). Other electronic 
orbitals do form bands which cross the Fermi level however, hence the metallic char- 
acter of the lanthanides. When pressure is applied, the f-electrons become increasingly 
itinerant. In fact, at some critical pressure, some rare-earth metals (mots notably Ce and 
Pr) undergo a sharp first-order transition which is accompanied by a discontinuous drop 
of the equilibrium unit-cell volume. Cerium is a particularly remarkable case, with a 
volume drop of as much as 15% and the same crystal symmetry (fee) in the low-volume 
(a) and high-volume (7) phase. In other cases, the transition corresponds to a change in 
crystal symmetry, from a lower symmetry phase at low pressure to a higher symmetry 
phase at high pressure. For a recent review on the volume-collapse transition of rare 
earth metals, see Ref. [15]. 

The equilibrium volume of actinide (5f) metals display behaviour which is interme- 
diate between transition metals and rare earths. From the beginning of the series (Th) 
until Plutonium (Pu), the volume has an approximately parabolic dependence on the 
filling of the f-shell, indicating delocalised 5f electrons. From Americium onwards, the 
volume has a much weaker dependence on the number of f-electrons, suggesting lo- 
calised behaviour. Interestingly, plutonium is right on the verge of this delocalisation to 
localisation transition. Not surprisingly then, plutonium is, among all actinide metals, 
the one which has the most complex phase diagram and which is also the most diffi- 
cult to describe using conventional electronic structure methods (see [16, 21] for recent 
reviews). This will be discussed further in the last section of these lectures. This very 
brief discussion of rare-earths and actinide compounds is meant to illustrate the need 
for methods able to deal simultaneously with the itinerant and localised character of 
electronic degrees of freedom. 



The physics of strong electronic correlations becomes even more apparent for f- 
electron materials which are compounds involving rare-earth (or actinide) ions and other 
atoms, such as e.g CeAls . A common aspect of such compounds is the formation of 
quasiparticle bands with extremely large effective masses (and hence large values of 
the low-temperature specific heat coefficient 7 = C/T), up to a thousand time the bare 
electron mass ! Hence the term "heavy-fermion" given to these compounds: for reviews, 
see e.g [22, 23]. The origin of these large effective masses is the weak hybridization be- 
tween the very localised f-orbitals and the rather broad conduction band associated with 
the metallic ion. At high temperature/energy, the f-electron have localised behaviour 
(yielding e.g local magnetic moments and a Curie law for the magnetic susceptibility). 
At low temperature/energy, the conduction electrons screen the local moments, leading 
to the formation of quasiparticle bands with mixed f- and conduction electron charac- 
ter (hence a large Fermi surface encompassing both f- and conduction electrons). The 
low-temperature susceptibility has a Pauli form and the low-energy physics is, apart 
from some specific compounds, well described by Fermi liquid theory. This screening 
process, the Kondo effect, is associated with a very low energy coherence scale, the 
(lattice [24]-) Kondo temperature, considerably renormalised as compared to the bare 
electronic energy scales. 
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FIGURE 3. Pressure volume data for the rare earths. Structures are identified, with "cmplx" signifiying 
a number of complex, low-symmetry structures. The volume collapse transitions are marked by the wide 
hatched lines for Ce, Pr, and Gd, while lines perpendicular to the curves denote the d-fcc to hP3 symmetry 
change in Nd and Sm. The curves are guides to the eye. Note that the data and curves have been shifted 
in volume by the numbers (in A^/atom) shown at the bottom of the figure. Figure and caption reproduced 
fromRef. [15]. 



The periodic Anderson model. The simplest model hamiltonian appropriate for f- 
electron materials is the Anderson lattice or periodic Anderson model. It retains the 

f-orbitals associated with the rare-earth or actinide atoms at each lattice site, as well as 
the relevant conduction electron degrees of freedom which hybridise with those orbitals. 
In the simplest form, the hamiltonian reads: 

HPAM = Lek4<7Ckc7 + £ {Vkcljmka+h.c) + £/ + L I E4(jm 

kcj kam Ram R \am j 

(V) 

Depending on the material considered, other terms may be necessary for increased 
realism, e.g an orbital dependent f-level hybridisation Vkm or interaction matrix 

^mm> ^ ^Vi^C'^ f-f hopping tff. 



2. DYNAMICAL MEAN-FIELD THEORY AT A GLANCE 

Dealing with strong electronic correlations is a notoriously difficult theoretical problem. 
From the physics point of view, the difficulties come mainly from the wide range of 
energy scales involved (from the bare electronic energies, on the scale of electron- 
Volts, to the low-energy physics on the scale of Kelvins) and from the many competing 
orderings and instabilities associated with small differences in energy. 
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FIGURE 4. Mean-field theory replaces a lattice model by a single site coupled to a self-consistent bath. 



It is the opinion of the author that, on top of the essential guidance from physical 
intuition and phenomenology, the development of quantitative techniques is essential in 
order to solve the key open questions in the field (and also in order to provide a deeper 
understanding of some "classic" problems, only partially understood to this day). 

In this section, we explain the basic principles of Dynamical Mean-Field Theory 
(DMFT). This approach has been developed over the last fifteen years and has led to 
some significant advances in our understanding of strong correlations. In this section, we 
explain the basic principles of this approach in a concise manner. The Hubbard model is 
taken as an example. For a much more detailed presentation, the reader is referred to the 
available review articles [3, 4]. 



2.1. The mean-field concept, from classical to quantum 

Mean-field theory approximates a lattice problem with many degrees of freedom by 
a single-site ejfective problem with less degrees of freedom. The underlying physical 
idea is that the dynamics at a given site can be thought of as the interaction of the local 
degrees of freedom at this site with an external bath created by all other degrees of 
freedom on other sites (Fig. 4). 



Classical mean-field theory. The simplest illustration of this idea is for the Ising 
model: 

// = -£7„5,5,-/i£5, (8) 

Let us focus on the thermal average of the magnetization on each lattice site: m, = (5,) . 
We consider an equivalent problem of independent spins: 

Heff = -l,hfs, (9) 
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in which the (Weiss) effective field is chosen in such a way that the value of m; is 
accurately reproduced. This requires: 

Phf^ = tanh-^mi (10) 

Let us consider, for definiteness, a ferromagnet with nearest-neighbour couplings 7,y = 
/ > 0. The mean-field theory approximation (first put forward by Pierre Weiss, under 
the name of "molecular field theory") is that h^/^ can be approximated by the thermal 
average of the local field seen by the spin at site i, namely: 

hf^c:th + Y,Jijmj = h + zJm (11) 
j 

where z is the connectivity of the lattice, and translation invariance has been used (Jij = J 
for n.n sites, mi = m). This leads to a self-consistent equation for the magnetization: 

m ^ tanh{l5h + zpJm) (12) 

We emphasize that replacing the problem of interacting spins by a problem of non- 
interacting ones in a effective bath is not an approximation, as long as we use this equiv- 
alent model for the only purpose of calculating the local magnetizations. The approxi- 
mation is made when relating the Weiss field to the degrees of freedom on neighbour- 
ing sites, i.e in the self-consistency condition (11). We shall elaborate further on this 
point of view in the next section, where exact energy functional will be discussed. The 
mean-field approximation becomes exact in the limit where the connectivity z of the 
lattice becomes large. It is quite intuitive indeed that the neighbors of a given site can 
be treated globally as an external bath when their number becomes large, and that the 
spatial fluctuations of the local field become negligible. 

Generalisation to the quantum case: dynamical mean-field theory. This construc- 
tion can be extended to quantum many-body systems. Key steps leading to this quan- 
tum generalisation where: the introduction of the limit of large lattice coordination for 
interacting fermion models by Metzner and VoUhardt [25] and the mapping onto a self- 
consistent quantum impurity by Georges and Kotliar [26], which established the DMFT 
framework^. 

I explain here the DMFT construction on the simplest example of the Hubbard 
modeP: 

H=-Y, ^ij j<^ + ^% "n"U + £0 1^ «!(T (1 3) 

ij,a i ia 

As explained above, it describes a collection of single- orbital "atoms" placed at the 
nodes R, of a periodic lattice. The orbitals overlap from site to site, so that the fermions 



See also the later work in Ref. [27], and Ref. [3] for an extensive list of references. 
^ The energy eo of the single-electron atomic level has been introduced in this section for the sake of 
pedagogy. Naturally, in the single band case, everything depends only on the energy Co — M with respect 
to the global chemical potential so that one can set eo = 
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can hop with an amplitude tij. In the absence of hopping, each "atom" has 4 eigenstates: 
|0), I t), I i) and ll) with energies 0, Sq and U + 2eo, respectively. 

The key quantity on which DMFT focuses is the local Green's function at a given 
lattice site: 

Gfi{T-T')^-{Tcia{T)cl{T')) (14) 

In classical mean-field theory, the local magnetization m, is represented as that of a single 
spin on site / coupled to an effective Weiss field. In a completely analogous manner, we 
shall introduce a representation of the local Green's function as that of a single atom 
coupled to an effective bath. This can be described by the hamiltonian of an Anderson 
impurity model ^: 

HaIM = Hatom + ^tef /i + ^coupling (15) 

in which: 

Hatom = U n'^n''^ + (Co - /i) + np 

Hcoupling = L/C7 {Aa^c + C^aia) (16) 

In these expressions, a set of non-interacting fermions (described by the aj's) have been 
introduced, which are the degrees of freedom of the effective bath acting on site R,. 
The £/ and V/'s are parameters which should be chosen in such a way that the c-orbital 
(i.e impurity) Green's function of (16) coincides with the local Green's function of the 
lattice Hubbard model under consideration. In fact, these parameters enter only through 
the hybridisation function: 

This is easily seen when the effective on-site problem is recast in a form which does 
not explicitly involves the effective bath degrees of freedom. However, this requires the 
use of an effective action functional integral formalism rather than a simple hamiltonian 
formalism. Integrating out the bath degrees of freedom one obtains the effective action 
for the impurity orbital only under the form: 

Seff = - f dT f dT'Y,ct{T)%\T-T')Ca{T') + U f dTn^{T)n^{T) (18) 
J JO Q JO 



in which: 



%\i(On) =i(On + ll-£o-Hi(On) (19) 



Strictly speaking, we have a collection of independent impurity models, one at each lattice site. In this 
section, for simplicity, we assume a phase with translation invariance and focus on a particular site of the 
lattice (we therefore drop the site index for the impurity orbital c^). We also assume a paramagnetic phase. 
The formalism easily generalizes to phases with long-range order (i.e translational and/or spin-symmetry 
breaking) [3] 
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This local action represents the effective dynamics of the local site under consideration: 
a fermion is created on this site at time T (coming from the "external bath", i.e from the 
other sites of the lattice) and is destroyed at time x' (going back to the bath). Whenever 
two fermions (with opposite spins) are present at the same time, an energy cost U is 
included. Hence this effective action describes the fluctuations between the 4 atomic 
states |0) , I I) , I i) , Ti) induced by the coupling to the bath. We can interpret %(t — i') 
as the quantum generalisation of the Weiss effective field in the classical case. The main 
difference with the classical case is that this "dynamical mean-field" is a. function of 
energy (or time) instead of a single number. This is required in order to take full account 
of local quantum fluctuations, which is the main purpose of DMFT. % also plays the 
role of a bare Green's function for the effective action Seff, but it should not be confused 
with the non-interacting {U — 0) local Green's function of the original lattice model. 

At this point, we have introduced the quantum generalisation of the Weiss effective 
field and have represented the local Green's function G„- as that of a single atom coupled 
to an effective bath. This can be viewed as an exact representation, as further detailed in 
Sec. 3. We now have to generalise to the quantum case the mean-field approximation 
relating the Weiss function to Gu (in the classical case, this is the self-consistency 
relation (12)). The simplest manner in which this can be explained - but perhaps not 
the more illuminating one conceptually (see Sec. 3 and [3, 28])- is to observe that, in 
the effective impurity model (18), we can define a local self-energy from the interacting 
Green's function G(t — t') = — < Tc{i)c'^{i') >Seff the Weiss dynamical mean- 
field as: 

^imp{i(On) = %^{i(On) - G-\iCO„) 

= i(On + H-eo-^{io)n)-G-^{io)n) (20) 

Let us, on the other hand, consider the self-energy of the original lattice model, defined 
as usual from the full Green's function Gij{T — t') = — < Tci^(j{t)c'j^{T') > by: 

G(k,ia)„) = i — (21) 

J(0„ + ;U - So - Ek - S(k, ««„) 

in which Ck is the Fourier transform of the hopping integral, i.e the dispersion relation 
of the non-interacting tight-binding band: 

ek = (22) 

j 

We then make the approximation that the lattice self-energy coincides with the impurity 
self-energy. In real-space, this means that we neglect all non-local components of L/^ 
and approximate the on-site one by 'Limp- 

^ii — ^imp ■) ^i^j — ^ (23) 

We immediately see that this is a consistent approximation only provided it leads to a 
unique determination of the local (on-site) Green's function, which by construction is 
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TABLE 1. Correspondance between the mean-field theory of a classical system and the dynamical 
mean-field theory of a quantum system. 



Quantum Case 


Classical Case 






H = — JijSiS j — hY,iSi 


Hamiltonian 


Gu{i(On) = -< cf{i(0„)ci{i(0„) > 


Mi =< Si > 


Local Observable 


^eff = Hatom + L/a £ia^a^la+ 
+ I/aV/(fl,+ C<;j + /l.c) 


^eff = -heffS 


Effective single-site 
Hamiltonian 


% ^ {i(0„) = i(On + A{i(0„) 


heff 


Weiss function/Weiss field 




Self-consistency relation 



the impuiity-model Green's function. Summing (21) over k in order to obtain the on- 
site component Gu of the the lattice Green's function, and using (20), we arrive at the 
self-consistency condition^: 

?A(,»,) + G(k)-'-% = 
Defining the non-interacting density of states: 

Z)(e) = £5(e-ek) (25) 

k 



this can also be written as: 



de , ^J:^^ , , = G(i(On) (26) 

This self-consistency condition relates, for each frequency, the dynamical mean-field 
A(/Ct)„) and the local Green's function G{iC0n). Furthermore, G(j(0„) is the interacting 
Green's function of the effective impurity model (16) -or (18)-. Therefore, we have a 
closed set of equations that fully determine in principle the two functions A, G (or %,G)). 
In practice, one will use an iterative procedure, as represented on Fig. 5. In many cases, 
this iterative procedure converges to a unique solution independently of the initial choice 
of A(ztt)„) . In some cases however, more than one stable solution can be found (e.g close 
to the Mott transition, see section below). The close analogy between the classical mean- 
field construction and its quantum (dynamical mean-field) counterpart is summarized in 
Table 1. 



^ Throughout these notes, the sums over momentum are normalized by the volume of the Brillouin zone, 
i.elkl = l 
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EFFECTIVE LOCAL IMPURITY PROBLEM 




FIGURE 5. The DMFT iterative loop. The following procedure is generally used in practice: starting 
from an initial guess for the impurity Green's function G,„,p is calculated by using an appropriate 
solver for the impurity model (top arrow). The impurity self-energy is also calculated from E„„p = 
^Q^{iO)„) — G^jlp{i(0„). This is used in order to obtain the on-site Green's function of the lattice model by 
performing a k-summation (or integration over the free d.o.s): G/oc = Lkl'^^^n + M ~ £k^'^iinp{i(On)]^^ ■ An 
updated Weiss function is then obtained as %^„g^, = Gy^^^ + T.i,„p, which is injected again into the impurity 
solver (bottom arrow). The procedure is iterated until convergence is reached. 



2.2. Limits in which DMFT becomes exact 

Two simple limits: non-interacting band and isolated atoms. It is instructive to check 
that the DMFT equations yield the exact answer in two simple limits: 

• In the non-interacting limit U = 0, solving (18) yields G{iCOn) = '^o{ico„) and 
'^imp = 0. Hence, from (24), G{i(On) = Ek 1/ (^^n + M — £o — Ck) reduces to the free 
on-site Green's function. DMFT is trivially exact in this limit since the self-energy 
is not only k-independent but vanishes altogether. 

• In the atomic limit tjj = 0, one just has a collection of independent atoms on each 
site and = 0. Then (24) implies A(?'(U„) = 0: as expected, the dynamical mean- 
field vanishes since the atoms are isolated. Accordingly, the self-energy only has 
on-site components, and hence DMFT is again exact in this limit. The Weiss field 
reads — i(0,i + jU — Cq, which means that the action Sg/f simply corresponds to 
the quantization of the atomic hamiltonian Hatom- This yields: 

r^l-f^ \ \—n/2 . n/2 

'^KlOJnjatom - + /(B„ +/!-(/ 

H^UJnjatom- 2 ^ ico„+ll-{l-n/2)U ' > 

with = ju -£o and n/2 = (^/^M +g/3(2M-f/))/(i +2e^^ 
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Hence, the dynamical mean-field approximation is exact in the two limits of the non- 
interacting band and of isolated atoms, and provides an interpolation in between. This 
interpolative aspect is a key to the success of this approach in the intermediate coupling 
regime. 

Infinite coordination. The dynamical mean-field approximation becomes exact in 
the limit where the connectivity z of the lattice is taken to infinity. This is also true 
of the mean-field approximation in classical statistical mechanics. In that case, the 
exchange coupling between nearest-neighbour sites must be scaled as: Jij = J /z (for 
Jifs of uniform sign), so that the Weiss mean-field hgff in (11) remains of order one. 
This also insures that the entropy and internal energy per site remain finite and hence 
preserves the competition which is essential to the physics of magnetic ordering. In the 
case of itinerant quantum systems [25], a similar scaling must be made on the hopping 
term in order to maintain the balance between the kinetic and interaction energy. The 
nearest-neighbour hopping amplitude must be scaled as: tij = tj^. This insures that 
the non-interacting d.o.s T>{e) = Lk ^(^ ~ ^k) has a non-trivial limit as z — > oo. Note that 
it also insures that the superexchange Jij oc tfj/U scales as 1/z, so that magnetic ordering 
is preserved with transition temperatures of order unity. In practice, two lattices are often 
considered in thez = °° limit: 

• The d-dimensional cubic lattice with z = 2d ^ o° and Ck = —'2tY,p=i cos(/:p) / y^. 
In this case the non-interacting d.o.s becomes a Gaussian: D{e) = 7^^^^P~(^) 

• The Bethe lattice (Cayley tree) with coordination z — > oo and nearest-neighbor hop- 
ping tij — t/y/z. This corresponds to a semicircular d.o.s: D{e) = ^y^l — (e/D)^ 
with a half-bandwidth D = 2t. In this case, the self -consistency condition (24) 
can be inverted explicitly in order to relate the dynamical mean-field to the local 
Green's function as: A(z(0„) = ?^G(ztt)„). 

Apart from the intrinsic interest of solving strongly correlated fermion models in the 
limit of infinite coordination, the fact that the DMFT equations become exact in this limit 
is important since it guarantees, for example, that exact constraints (such as causality 
of the self-energy, positivity of the spectral functions, sum rules such as the Luttinger 
theorem or the f-sum rule) are preserved by the DMFT approximation. 



2.3. Important topics not reviewed here 

There are several important topics related to the DMFT framework, which I have not 
included in these lecture notes. Some of them were covered in the lectures, but extensive 
review articles are available in which these topics are at least partially described. 

This is a brief list of such topics: 

DMFT for ordered phases. The DMFT equations can easily be extended to study 
phases with long-range order, calculate critical temperatures for ordering as well as 
phase diagrams, see e.g [3]. 
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Response and correlation functions in DMFT. Response and correlation functions 
can be expressed in terms of the lattice Green's functions, and of the impurity model 
vertex functions, see e.g [3, 4]. Note that momentum-dependence enters, through the 
lattice Green's function. 

Physics of the Anderson impurity model. Understanding the various possible fixed 
points of quantum impurity models is important for gaining physical intuition when 
solving lattice models within DMFT. See Ref. [23] for a review and references on the 
Anderson impurity model. It is important to keep in mind that, in contrast to the common 
situation in the physics of magnetic impurities or mesoscopics, the effective conduction 
electron bath in the DMFT context has significant energy-dependence. Also, the self- 
consistency condition can drive the effective impurity model from one kind of low- 
energy behaviour to another, depending on the range of parameters (e.g close to the 
Mott transition, see Sec. 4). 

Impurity solvers. Using reliable methods for calculating the impurity Green's func- 
tion and self-energy is a key step in solving the DMFT equations. A large numbers 
of "impurity solvers" have been implemented in the DMFT context^, including: the 
quantum Monte Carlo (QMC) method [27] (see also [29, 30]), based on the Hirsch- 
Fye algorithm [31], adaptative exact diagonalisation or projective schemes (see [3] for a 
review and references), the Wilson numerical renormalisation group (NRG, see e.g [32] 
and references therein). Approximation schemes have also proven useful, when used 
in appropriate regimes, such as the "iterated perturbation theory" approximation (IPT, 
[26, 33]), the non-crossing approximation (NCA, see [4] for references) and various 
extensions [34], as well as schemes interpolating between high and low energies [35]. 

Beyond DMFT. DMFT does capture ordered phases, but does not take into account 
the coupling of short-range spatial correlations (let alone long- wavelength) to quasiparti- 

cle properties, in the absence of ordering. This is a key aspect of some strongly correlated 
materials (e.g cuprates, see the concluding section of these lectures), which requires an 
extension of the DMFT formalism. Two kinds of extensions have been explored: 

• k-dependence of the self-energy can be reintroduced by considering cluster exten- 
sions of DMFT, i.e a small cluster of sites (or coupled atoms) into a self-consistent 
bath. Various embedding schemes have been discussed [3, 36, 37, 38, 39, 40, 41] 
and I will not attempt a review of this very interesting line of research here. One of 
the key questions is whether such schemes can account for a strong variation of the 
quasiparticle properties (e.g the coherence scale) along the Fermi surface. 

• Extended DMFT (E-DMFT [42, 43, 44, 45]) focuses on two-particle local observ- 
ables, such as the local spin or charge correlation functions, in addition to the lo- 
cal Green's function of usual DMFT. For applications to electronic structure, see 
Sec. 5.5. 



Some early versions of numerical codes are available at: http://www.lps.ens.fr/~krauth 
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3. FUNCTIONALS, LOCAL OBSERVABLES, AND 
INTERACTING SYSTEMS 



In this section^, I would like to discuss a theoretical framework which applies quite 

generally to interacting systems. This framework reveals common concepts underlying 
different theories such as: the Weiss mean-field theory (MFT) of a classical magnet, 
the density functional theory (DFT) of the inhomogeneous electron gas in solids, and 
the dynamical mean-field theory (DMFT) of strongly correlated electron systems. The 
idea which is common to these diverse theories is the construction of a functional of 
some local quantity (effective action) by the Legendre transform method. Though exact 
in principle, it requires in practice that the exact functional is approximated in some 
manner. This method has a wide range of applicability in statistical mechanics, many- 
body physics and field-theory [46]. The discussion will be (hopefully) pedagogical, and 
for this reason I will begin with the example of a classical magnet. For a somewhat more 
detailed presentation, see Ref. [28]. 

There are common concepts underlying all these constructions (cf. Table), as will 
become clear below, namely: 

• i) These theories focus on a specific local quantity: the local magnetization in MFT, 
the local electronic density in DFT, the local Green's function (or spectral density) 
in DMFT. 

• ii) The original system of interest is replaced by an equivalent system, which is used 
to provide a representation of the selected quantity: a single spin in an effective field 
for a classical magnet, free electrons in an effective one-body potential in DFT, a 
single impurity Anderson model within DMFT. The effective parameters entering 
this equivalent problem define generalized Weiss fields (the Kohn-Sham potential 
in DFT, the effective hybridization within DMFT), which are self-consistently 
adjusted. I note that the associated equivalent system can be a non-interacting (one- 
body) problem, as in MFT and DFT, or a fully interacting many-body problem 
(albeit simpler than the original system) such as in DMFT and its extensions. 

• iii) In order to pave the way between the real problem of interest and the equivalent 
model, the method of coupling constant integration will prove to be very useful 
in constructing (formally) the desired functional using the Legendre transform 
method. The coupling constant can be either the coefficient of the interacting part of 
the hamiltonian (which leads to a non-interacting equivalent problem, as in DFT), 
or in front of the non-local part of the hamiltonian (which leads in general to a 
local, but interacting, equivalent problem such as in DMFT). 

Some issues and questions are associated with each of these points: 

• i) While the theory and associated functional primarily aims at calculating the se- 
lected local quantity, it always come with the possibility of determining some more 
general object. For example, classical MFT aims primarily at calculating the local 
magnetization, but it can be used to derive the Ornstein-Zemike expression of the 



^ This section is based in part on Ref. [28] 
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TABLE 2. Comparison of theories based on functionals of a local observable 



Theory 


MFT 


DFT 


DMFT 


Quantity 


Local magnetization m, 


Local density n{x) 


Local GFG;,(co) 


Equivalent 
system 


Spin in 
effective field 


Electrons in 
effective potential 


Quantum 

impurity model 


Generalised 
Weiss field 


Effective 
local field 


Kohn-Sham 
potential 


Effective 
hybridisation 



correlation function between different sites. Similarly, DFT aims at the local den- 
sity, but Kohn-Sham orbitals can be interpreted (without a firm formal justification) 
as one-electron excitations. DMFT produces a local self-energy which one may in- 
terpret as the lattice self -energy from which the full k-dependent Green's function 
can be reconstructed. In each of these cases, the precise status and interpretation of 
these additional quantities can be questioned. 

• ii) I emphasize that the choice of an equivalent representation of the local quan- 
tity has nothing to do with subsequent approximations made on the functional. The 
proposed equivalent system is in fact an exact representation of the problem under 
consideration (for the sake of calculating the selected local quantity). It does raise 
a representability issue, however: is it always possible to find values of the gen- 
eralised Weiss field which will lead to a specified form of the local quantity, and 
in particular to the exact form associated with the specific system of interest? For 
example: given the local electronic density n{x) of a specific solid, can one always 
find a Kohn-Sham effective potential such that the one-electron local density ob- 
tained by solving the Schrodinger equation in that potential coincides with n{x) ? 
Or, in the context of DMFT: given the local Green's function of a specific model, 
can one find a hybridisation function such that it can be viewed as the local Green's 
function of the specified impurity problem ? 

• iii) There is also a stability issue of the exact functional: is the equilibrium value of 
the local quantity a minimum ? More precisely, one would like to show that negative 
eigenvalues of the stability matrix correspond to true physical instabilities of the 
system. I will not seriously investigate this issue in this lecture (for a discussion 
within DMFT, where it is still quite open, see [47]). 



3.1. The example of a classical magnet 

For the sake of pedagogy, I will consider in this section the simplest example on which 
the above ideas can be made concrete: that of a classical Ising magnet with hamiltonian 

H=-Y,JijSiSj (28) 

ij 

Construction of the effective action. We want to construct a functional r[m,] of a 
preassigned set of local magnetizations m,, such that minimizing this functional yields 
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the equilibrium state of the system. This functional is of course the Legendre transform 
of the free-energy with respect to a set of local magnetic fields. To make contact with 
the field-theory literature, I note that j8 F is generally called the effective action in this 
context. I will give a formal construction of this functional, following a method due to 
Plefka [48] and Yedidia and myself [49]. Let us introduce a varying coupling constant 
a e [0, 1], and define: 

Ha = aH = Y^aJijSiSj (29) 

ij 

Introducing local Lagrange multipliers A,-, we consider the functional: 

Q.[miMa\ = -^lnTre-/^««+^^'^'(^'-'«') =F[A,]+J2A;m; (30) 

Requesting stationarity of this functional with respect to the A,'s amounts to impose 
that, /or all values of a, (Sj) coincides with the preassigned local magnetization m,. The 
equations m,- = (S,) which expresses the magnetization as a function of the sources A/ 
can then be inverted to yield the A;'s as functions of the m/s and of a: 

{Si)x^a = m^k^h[mi;a\ (31) 

(The average (■ ■ ■)x^a ™^ this equation is with respect to the Boltzmann weight appearing 
in the above definition of fl, including A/s and a). The Lagrange parameters can then 
be substituted into Q. to obtain the a-dependent Legendre transformed functional: 

Ta [mi] = ^ [mi , A, [mj ,a]]=F [A, [m] ] + A/ [m] (32) 

Of course, the functional we are really interested in is that of the original system with 
a = 1, namely: 

r[mi]=ra=i[mi] (33) 

Let us first look at the non-interacting limit a = for which the explicit expression of 
Q. is easily obtained as: 

= £ (^-^IncoshjSAi + m/A;^ (34) 

Varying in the A's yields: 

tanhj8A/"=°^=m/ (35) 



and finally: 



a 



1^/1+m,, l+mi , l-m,^^l-m, ^^^^ 



1 \- / 1+m; l+nii 



The a = theory defines the equivalent problem that we want to use in order to deal 
with the original system. Here, it is just a theory of independent spins in a local effective 
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field. The expression (36) is simply the entropy term corresponding to independent Ising 
spins for a given values of the local magnetizations. 

The value taken by the Lagrange multiplier in the equivalent system, A, """^ (denoted 
Xf in the following), must be interpreted as the Weiss effective field. We note that, in this 
simple example, there is an explicit and very simple relation (35) between the Weiss field 
and m,-, so that one can work equivalently in terms of either quantities. Also, because of 
the simple form of (35), representability is trivially satisfied: given the actual values of 
the magnetizations m,'s (e [—1,1]) at equilibrium for the model under consideration, 

one can always represent them by the Weiss fields Ph^/^ = arctanhm;. 

To proceed with the construction of F, we use a coupling constant integration and 
write: 

TK; ce = 1] = ToK] + / da—^[mi] (37) 

Jo da 

It is immediate that, because of the constraint {{Si — mi)) = 0: 

-J^ = {H)a,X[a] = - '^Jij{SiSj)a,X[a,m] (38) 

In this expression, the correlation must be viewed as a functional of the local magne- 
tizations (thanks to the inversion formula (31)). Introducing the connected correlation 
function: 

g^ij[{mk};a] = {{Si-mi){Sj-mj))a^x[a,m] (39) 

we obtain: 

-7^ = - Y^Jijtnmj - g'ij [{mk} ; «] (40) 
^" ij ij 

So that finally, one obtains the formal expression for r[mj] = ra=i [mi] : 

r[m/] = ro[m/] - yjijtnimj - YJij I dagUmu; «] = Tq + Emf + ^corr (41) 
ij ij 

In this expression, denotes the connected correlation function for a given value of the 
coupling constant, expressed as a fiinctional of the local magnetisations. 
Hence, the exact functional T appears as a sum of three contributions: 

• The part associated with the equivalent system (corresponding here to the entropy 
of constrained but otherwise free spins) 

• The mean-field energy Y^ij^ij^i^j 

• A contribution from correlations which contains all corrections beyond mean-field 

As explained in the next section, there is a direct analogy between this and the various 
contributions to the density functional within DFT (kinetic energy, Hartree energy and 
exchange-correlation) . 

I note in passing that one can derive a closed equation for the exact functional, which 
reads (see [28] for a derivation): 
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1 

Talmi] = Tolmi] - aJ^Jijniimj - - J^Jij da' 



5m]^5mi 



(42) 

-I I] 



This equation fully determines in principle the effective action functional. However, in 
order to use it in practice, one generally has to start from a limit in which the functional 
is known explicitly, and expand around that limit. For example, an expansion around the 
high-temperature limit yields systematic corrections to mean-field theory [49, 28]. This 
equation is closely related [28] to the Wilson-Polchinsky equation [50] for the effective 
action (after a Legendre transformation: see also [51]), which can be taken as a starting 
point for a renormalisation group analysis by starting from the local limit and expanding 
in the "locality" (see e.g [52, 51]. 

Equilibrium condition and stability. The physical values of the magnetisations at 
equilibrium are obtained by minimising F, which yields: 

mr = tanh(^^£7,m;.-^^j (43) 

and the Weiss field takes the following value: 

/.f^^^(AP)* = p,m*-^r (44) 

This equation is a self-consistency condition which determines the Weiss field in terms 
of the local magnetizations on all other sites. Its physical interpretation is clear: h^^-^ is 
the true (average) local field seen by site i. It is equal to the sum of two terms: one in 
which all spins are treated as independent, and a correction due to correlations. 

The stability of the functional around equilibrium is controlled by the fluctuation 
matrix: 

^!l_ = M_y^..+i!E^ (45) 

dniidnij 8m j dmidmj 

At equilibrium, this is nothing else than the inverse of the susceptibility (or correlation 
function) matrix: 

iX ')ij = iXo')ij-Jij + T-^. (46) 



5mi5mj o v y dmiSmj 

with: 

(^«"')« = mT^''^ '''' 

Hence, our functional does satisfy a stability criterion as defined in the introduction: a 
negative eigenvalue of this matrix (i.e of x{^)) would correspond to a physical instability 
of the system. Note that at the simple mean-fied level, we recover the RPA formula for 
the susceptibility: {x~^)ij = iXo^)ij-Jij- 
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Mean-field approximation and beyond. Obviously, this construction of the exact 
Legendre transformed free energy, and the exact equilibrium condition (43) has formal 
value, but concrete applications require some further approximations to be made on 
the correlation term Tcorr- The simplest such approximation is just to neglect Tcorr 
altogether. This is the familiar Weiss mean-field theory: 

1 /1+m/ l+m; 1-m; 1-mA ^ 
^MFT = L I + J ~ l^Jijmmj (48) 



For a ferromagnet (uniform positive 7//s), this approximation becomes exact in the limit 
of infinite coordination of the lattice. 

The formal construction above is a useful guideline when trying to improve on the 
mean-field approximation. I emphasize that, within the present approach, it is the self- 
consistency condition (44) (relating the Weiss field to the environment) that needs to be 
corrected, while the equation m, = tanhj8/ij^^ is attached to our choice of equivalent 
system and will be always valid. For example, in [48, 49] it was shown how to con- 
struct Tcorr by a systematic high-temperature expansion in j8. This expansion can be 
conveniently generated by iterating the exact equation (42). It can also be turned into an 
expansion around the limit of infinite coordination [49]. The first contribution to Tcorr 
in this expansion appears at order j8 (or a^) and reads: 

TSrr = ~l,jfj{l-mj)il-mj) (49) 

^ ij 

This is a rather famous correction to mean-field theory, known as the "Onsager reaction 
term". For spin glass models (Ji/s of random sign), it is crucial to include this term 
even in the large connectivity limit. The corresponding equations for the equilibrium 
magnetizations are those derived by Thouless, Anderson and Palmer [53]. 



3.2. Density functional theory 

In this section, I explain how density-functional theory ^ (DFT) [55, 56] can be derived 
along very similar lines. This section borrows from the work of Fukuda et al. [57, 46] 
and of Valiev and Fernando [58]. For a recent pedagogical review emphasizing this point 
of view, see [59]. For detailed reviews of the DFT formalism, see e.g [60, 61]. 

Let us consider the inhomogeneous electron gas of a solid, with hamiltonian: 

^ = - E ^ V2 + l^vin) + ^ £ f/(r,- r,-) (50) 

i i i^j 

in which v{x) is the external potential due to the nuclei and U{x— x') (= — |) is 
the electron-electron interaction. (I use conventions in which h — m— 1). Let us write 



° I actually consider the finite-temperature extension of DFT [54] 
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this hamiltonian in second- quantized form, and again introduce a coupling-constant 
parameter a (the physical case is a — 1): 

Ha = ~Jdx\if^V^\}f+ J dxv{x)n{x) + ^ J dxdx! n{x)U{x-x')h{x') (51) 

We want to construct the free energy functional of the system while constraining the 
average density to be equal to some specified function n{x). In complete analogy with 
the previous section, we introduce a Lagrange multiplier function A (x), and consider ^: 

t2„[n(;c),A(;c)] = -^lnTrexp(^-j8//a + /3 j dxX{x){n{x)-h{x))^ (52) 

A functional of both n{x) and As before, stationarity in A insures that: 

(«W)A,a = n{x) X{x)^ K[n{x)\ (53) 

This will be used to eliminate A (x) in terms of n{x) and construct the functional of n{x) 
only: 

Ta[n{x)]=^aHx),laHx)]] (54) 



3.2.1. Equivalent system: non-interacting electrons in an effective 

potential 

Again, I first look at the non-interacting case a = 0. Then we have to solve a one- 
particle problem in an x-dependent external potential. This yields: 

^2o['^H:■^M] = —tvln[iO}n — i — v — X] — J dxX{x)n{x) (55) 

In this equation, tr denotes the trace over the degrees of freedom of a single electron, i(On 
is the usual Matsubara frequency, and f = — v, A are the one-body operators cor- 
responding to the kinetic energy, external potential and X (x) respectively. The identity 
In det = tr In has been used. 

Minimisation with respect to A (x) yields the following relation between and n{x) : 

P „ lO^-f-v-Ao 

This defines the functional Aq albeit in a somewhat implicit manner. This is 

directly analogous to Eq.(35) defining the Weiss field in the Ising case (but in that case. 



^ Note that I chose in this expression a different sign convention for A than in the previous section, and 
also that Tr denotes the full many-body trace over all A'-electrons degrees of freedom. 
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this equation was easily invertible). If we want to be more explicit, what we have to do 
is solve the one-particle Schrodinger equation: 



■^A + VKsix)^ Mx) = eMx) (57) 

where the effective one-body potential (Kohn-Sham potential) is defined as: 

VKsix)=v{x) + l\x) (58) 
It is convenient to construct the associated resolvent: 

R(,y.jc^) = imn^ (59) 

, icOn-ei 

and the relation (56) now reads: 

Ll0/W|VFZ)(e/)=n(x) (60) 

in which ffo is the Fermi-Dirac distribution. 

This relation expresses the local density in an interacting many-particle system as 
that of a one-electron problem in an effective potential defined by (56). In so doing, 
the effective one-particle wave functions and energies (Kohn-Sham orbitals) have been 
introduced, whose relation to the original system (and in particular their interpretation 
as excitation energies) is far from obvious (see e.g [61]). There is, for example, no 
fundamental justification in identifying the resolvent (59) with the true one-electron 
Green's function of the interacting system. The issue of representability (i.e whether an 
effective potential can always be found given a density profile n{x)) is far from being as 
obvious as in the previous section, but has been established on a rigorous basis [62, 63]. 

To summarize, the non-interacting functional Fq [n(jc)] reads: 

ro[n(x)] = — trln[/(a^, — ? — V — Ao[n]] — J dxX^[x;n\n{x) (61) 
which can be rewritten as: 

ro[n(jc)] == -4£ln [l _ f dxvKs{x)n{x) + I dxv{x)n{x) (62) 

in which Xq and v^-^ are viewed as a functional of n{x), as detailed above. 
In the limit of zero temperature (j8 oo), this reads: 

ro[n(x), r = 0] = ^^l~ J dxvKsix)n{x) + J dxv{x)n{x) (63) 

in which the sum is over the N occupied Kohn-Sham states. We note that it contains extra 
terms beyond the ground-state energy of the KS equivalent system (see also Sec. 5.4). 
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We also note that Tq is not a very explicit functional of n{x). It is a somewhat more 
explicit functional of Ao(^) (or equivalently of the KS effective potential vks{x)) so that 
it is often more convenient to think in terms of this quantity directly. At any rate, in 
order to evaluate Tq for a specific density profile or effective potential one must solve 
the Schrodinger equation for KS orbitals and eigenenergies. This is a time-consuming 
task for realistic three-dimensional potentials and practical calculations would be greatly 
facilitated if a more explicit accurate expression for rrn(;c)] would be available 



3.2.2. The exchange-correlation functional 

We turn to the interacting theory, and use the coupling constant integration method 
(see [64] for its use in DFT): 

r\ AT 

T[n{x)] = r[n(jc); a = 0] + / da—^ (64) 

Jo da 



{U)x,a = y dxdx'U{x-x'){n{x)h{x'))^^a (65) 



Similarly as before: 

dTg 
da 

Separating again a Hartree (mean-field) term, we get: 

r[n{x)] = ro[n{x)] + EHanree[n{x)] +rxc[n{x)] (66) 

with: ^ 

Enartree [n{x)] = 2 / x)n{x)n{x') (67) 

and Txc is the correction-to mean field term (the exchange-correlation functional): 

rxc[n{x)] = ^Jdxd:>^U{x-x') J dag'a[n;x,x'] (68) 

In which: 

ga[n;x,x] = {{n{x) -n{x)){n{x) -n{x)))i^[n]^a (69) 

is the (connected) density-density correlation function, expressed as a functional of the 
local density, for a given value of the coupling a. 

It should be emphasized that the exchange-correlation functional T^c is independent of 
the specific form of the crystal potential v{x): it is a universal functional which depends 
only on the form of the inter-particle interaction U{x— x') ! To see this, we first observe 
that, because r[n(x)] is the Legendre transform of the free energy with respect to the 



see e.g the lecture notes by K.Burke: http://dft.rutgers.edu/k;ieron/beta/index.html 
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one-body potential, we can easily relate the functional in the presence of the crystal 
potential v(jc) to that of the homogeneous electron gas (i.e with v = 0): 



'[n{x)] = rHEG[n{x)]+ J dxv{x)n{x) 



(70) 



Since this relation is also obeyed for the non-interacting system (see Eq. (61)), and using 
r = Fq + F// + Txc, we see that the functional form of F_^c is independent of v(jc). It is 
the same for all solids, and also for the homogeneous electron gas. 

I finally note that an exact relation can again be derived for the density functional (or 
alternatively the exchange-correlation functional) by noting that: 



pgl,[n;x,x'] = 



5r« 



5n{x)5n{y) 

Inserting this relation into (66,68), one obtains: 



(71) 



1 /■ 

'^a[n]=ro[n]+aEH[n] + -JdxdxU{x-x') J da' 



dn(x)dn(y) 



(72) 



xx' 



in complete analogy with (42). For applications of this exact functional equation, see 
e.g [65, 66]. Analogies with the exact renormalization group approach (see previous 
section) might suggest further use of this relation in the DFT context. 



3.2.3. The Kohn-Sham equations 

Let us now look at the condition for equilibrium. We vary F[n(x)], and we note that, 
as before, the terms originating from the variation 5n{x) cancel because of the 
relation (56). We thus get: 

= -Ao(x)+ / dx'U{x-x')n{x') + ^^ (73) 
on[x) J on[x) 

so that the equilibrium density n*{x) is determined by: 

X\xy = [ ^ix'[/(x-x')n*(x') + |^|„=„* (74) 
J on[x) 

which equivalently specifies the KS potential at equilibrium as: 

^Ksi^) ^ ^(•^) + jdxU {x-x)n* (x) + ln=n* C^^) 

Equation (74) is the precise analog of Eq.(44) determining the Weiss field in the Ising 

case, and v^^ is the true effective potential seen by an electron at equilibrium, in a one- 
electron picture. Together with (57), it forms the fundamental (Kohn-Sham) equations 
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of the DFT approach. To summarize, the expression of the total energy {T = 0) reads: 

r[n{x),T = 0] = ^g/ - y dxvKs{x)n{x) + J dxv{x)n{x) + Txc [n{x)] (76) 

Concrete appUcations of the DFT formaUsm require an approximation to be made on 
the exchange-correlation term. The celebrated local density approximation (LDA) reads: 

r^c [n{x)] \lda = J dxn{x) e^J'^i^)] (77) 

in which e^f''^{n) is the exchange-correlation energy density of the homogeneous elec- 
tron gas, for an electron density n. Discussing the reasons for the successes of this ap- 
proximation (as well as its limitations) is quite beyond the scope of these lectures. The 
interested reader is referred e.g to [61, 59]. 

Finally, we observe that DFT satisfies the stability properties discussed in the intro- 
duction, since 5^T / 5n{x)5n{x') is the inverse of the density-density response function 
(q-dependent compressibility). A negative eigenvalue would correspond to a charge or- 
dering instability. 



3.3. Exact functional of the local Green's function, and the Dynamical 

Mean-Field Theory approximation 

In this section, I would like to explain how the concepts of the previous sections pro- 
vide a broader perspective on the dynamical mean field approach to strongly correlated 
fermion systems. In contrast to DFT which focuses on ground-state properties (or ther- 
modynamics), the goal of DMFT (see [3] for a review) is to address excited states by 
focusing on the local Green's function (or the local spectral density). Thus, it is natural 
to formulate this approach in terms of a functional of the local Green's function. This 
point of view has been recently emphasized by Chitra and Kotliar [67] and by the author 
in Ref. [28]. 

I describe below how such an exact functional can be formally constructed for a 
correlated electron model (irrespective, e.g of dimensionality), hence leading to a local 
Green's function (or local spectral density) functional theory. I will adopt a somewhat 
different viewpoint than in [67], by taking the atomic limit (instead of the non-interacting 
limit) as a reference system. This leads naturally to represent the exact local Green's 
function as that of a quantum impurity model, with a suitably chosen hybridisation 
function. There is no approximation involved in this mapping (only a representability 
assumption). This gives a general value to the impurity model mapping of Ref. [26]. 
Dynamical mean field theory as usually implemented can then be viewed as a subsequent 
approximation made on the non-local contributions to the exact functional (e.g. the 
kinetic energy). 

For the sake of simplicity, I will take the Hubbard model as an example throughout 
this section. The hamiltonian is decomposed as: 

Ha = UY, ~ " L ^'7 4a<^ jo (78) 
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I emphasize that the varying coupling constant a G [0, 1] has been introduced in front 
of the hopping term, which is the non-local term of this hamiltonian, and not in front 
of the interaction. When dealing with a more general hamiltonian, we would similarly 
decompose H = H^^^ + a^non-loc 



3.3.1. Representing the local Green's function by a quantum impurity 

model 

In order to constrain the local Green's function (c/(t)c^(t')) to take a specified value 
G(t — t'), we introduce conjugate sources (or Lagrange multipliers) A(t — x') and 
consider 

ri8 



^a[G{(0),^{(0)] = -^1"! DcDc+Qxp{J^ dT(£c+{-d, + ^)cia - Ha[c,c+]) + 

+ f f dTdT'Y,K^-^')[G{T-T')-cl{T)ci„{T')]} (79) 



Inverting the relation G = Ga [A] yields A = [G] , and a functional of the local Green's 
function is obtained as Ta[G] = Q.a [G,Aa[G]]. This is the Legendre transform of the 
free energy with respect to the local source A. 

I would like to emphasize that this construction is quite different from the Baym- 
Kadanoff formalism, which considers a functional of all the components of the lattice 
Green's function G,/, not only of its local part Gu. The Baym-Kadanoff approach also 
gives interesting insights into the DMFT construction [3, 47], and will be considered at 
a later stage in these lectures. 

Consider first the a = case, in which the hamiltonian is purely local (atomic limit). 
Then, we have to consider a local problem defined by the action: 



Simp= - [ dT I dTX4(^)[(-^T + M)5(T-T')-Ao(T-T')]cc,(T') 
JO JO Q 

+U / dTn^{T)ni{T) (80) 

J 



Hence, the local Green's function G{iC0n) is represented as that of a quantum impurity 
problem (an Anderson impurity problem in the context of the Hubbard model): 

G^Gi^p[Ao] (81) 

As before, Aq plays the role of a Weiss field (analogous to the effective field for a magnet, 
or to the KS effective potential in DFT). Formally, this Weiss field specifies [26] the 



In this section, I will divide the free energy functional by the number A^, of lattice sites (restricting 
myself for simplicity to an homogeneous system) 
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effective bare Green's function of the impurity action (80): 



% ^ (icOn) = i(On + H-Aq {icd„) (82) 

There are however two important new aspects here: 

• i) The Weiss function Aq is a dynamical (i.e frequency dependent) object. As a 
resuh the local equivalent problem (80) is not in Hamiltonian form but involves 
retardation 

• ii) The equivalent local problem is not a one-body problem, but involves local 
interactions. 

We note that, as in DFT, the explicit inversion of (81) is not possible in general. In 
practice, one needs a (numerical or approximate) technique to solve the quantum im- 
purity problem (an "impurity solver"), and one can use an iterative procedure. Starting 
from some initial condition for Aq (or ^0) , one computes the interacting Green's func- 
tion Gimp, and the associated self-energy Li^p = — ^7mp- ^^en updates ^0 as: 
imp + G ^] ^, where G is the specified value of the local Green's function. 

3. 3. 2. Exact Junctional of the local Green 's function 

We proceed with the construction of the exact functional of the local Green's function, 
by coupling constant integration (starting from the atomic limit). 

At a = (decoupled sites, or infinitely separated atoms), we have : 

[Ao, G] = Fi,np [Ao] - Tr (GAq) (83) 

where Funp is the free energy of the local quantum impurity model viewed as a functional 
of the hybridisation function. By formal inversion Aq = Ao[G] : 

ro[G] = Fimp [Ao[G]] - Tr (GAo[G] ) (84) 

We then observe that (since the a -derivatives of the Lagrange multipliers do not con- 
tribute because of the stationarity of 

= (^non-loc) (^5) 
which, for the Hubbard model, reduces to the kinetic energy: 

^ = (3^) = -^l.tij{ctcj)\G = Tr ^i:ekG«(k,/o^)|G (86) 

ij k 

In this expression, the lattice Green's function Ga (k, i(On) should be expressed, for a 
given a, as a functional of the local Green's function G. 



In this formula and everywhere below, Tr denotes j 2^„, with possibly a convergence factor e"*^o+ 
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This leads to the following formal expression of the exact functional r[G] = ra=i [G] : 
r[G] = Fi„p[Ao[G]]-TT (GAo[G]) + ^[G] (87) 
in which ^[G] is the kinetic energy functional (evaluated while keeping G„- = G fixed): 

^[G] = ['da^'£tij{c+cj)\G= daTv ^Y.^i,Ga{KiC0n)\G (88) 

•^0 jj Jo ^^s j{ 

The condition dT/ 5G = determines the actual value of the local Green's function at 
equilibrium as (using STq/ 5G = — Aq): 

^i«('-'*)i = S '''' 

We recall that the generalized Weiss function (hybridization) and G are, by construction, 
related by (81): 

G = Gimpi^o] (90) 

Equations (89,90) (together with the definition of the impurity model, Eq. 80)) are 
the key equations of dynamical mean-field theory, viewed as an exact approach. The 
cornerstone of this approach [26] is that, in order to obtain the local Green's function, 
one has to solve an impurity model (80), submitted to the self-consistency condition (89) 
relating the hybridization function Aq to G{i(0„) itself. I emphasize that, since r[G] is 
an exact functional, this construction is completely general: it is valid for the Hubbard 
model in arbitrary dimensions and on an arbitrary lattice. 

Naturally, using it in practice requires a concrete approximation to the kinetic energy 
functional ^[G] (similarly, the DPT framework is only practical once an approximation 
to Txc is used, for example the LDA). The DMFT approximation usually employed is 
described below. In fact, it might be useful to employ a different terminology and call 
"local spectral density functional theory" (or "local impurity functional theory") the 
exact framework, and DMFT the subsequent approximation commonly made in ^[G] . 



3.3.3. A simple case: the infinite connectivity Bethe lattice 

It is straightforward to see that the formal expression for the kinetic energy functional 
^[G] simplifies into a simple closed expression for the Bethe lattice with connectivity 
z, in the limit z — > oo. In fact, a closed form can be given on an arbitrary lattice in the 
limit of large dimensions, but this is a bit more tedious and we postpone it to the next 
section. 

In the limit of large connectivity, the hopping must be scaled as: tij = tj ^ [25]. 
Expanding the kinetic energy functional in (87) in powers of oc, one sees that only the 
term of order a remains in the z = °° limit thanks to the tree-like geometry, namely: 

fi 

ocY,tijtki{4cj4ci)a=o = aY^tf/Tv g2 = a{zNs) -Tr G^ (91) 

ijkl ij ^ 
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So that, integrating over a, one obtains ^[G] = t^Tv and finally: 

rBethe,z=4G] = Fi„,p [Ao[G]] - Tr (GAo[G]) + -Tr (92) 

This functional is similar (although different in details) to the one recently used by 
Kotliar [68] in a Landau analysis of the Mott transition within DMFT. 

The self-consistency condition (89) that finally determines both the local Green's 
function and the Weiss field (through an iterative solution of the impurity model) thus 
reads in this case: 

Ao[G,icOn]=t^GiicOn) (93) 



3.3.4. DMFT as an approximation to the kinetic energy Junctional. 

Now, I will show that the usual form of DMFT [3] (for a general non-interacting dis- 
persion e^) corresponds to a very simple approximation of the kinetic energy term -^[G] 
in the exact functional r[G]. Consider the one-particle Green's function Ga(k,/tt)„) as- 
sociated with the action (79) of the Hubbard model, in the presence of the source term 
Aa and for an arbitrary coupling constant. We can define a self-energy associated with 
this Green's function: 

Ga (k, ico„) = (94) 

i(On + H- Aa [i(On] - ttCk - I^a [K iOh] 

The self-energy Ha is in general a k-dependent object, except obviously for a = in 
which all sites are decoupled into independent impurity models. The DMFT approxima- 
tion consists in replacing La for arbitrary a by the impurity model self-energy Eq (hence 
depending only on frequency), at least for the purpose of calculating £^[G]. Hence: 

Ga (k, ia)n) \dmft = ■ ^ r. ^ (95) 

iO}„ + li-Aa[iO}„;G]-aey,-La=o[iO}„,G] 



With: 



2:a=o[G;i(oJ =^0-1 -G-i = io}n + li-Ao[i(On.G]^G-^ (96) 



Summing over k, one then expresses the local Green's function in terms of the hybridi- 
sation as: 

G(i.,)=/*^ = i6(i) (97) 

With ^ = i(On + H — Aa —Lq = Aq — Aa + G~^ . In this expression, D(e) = 7^Lk^(£~£k) 
is the non-interacting density of states, and D{z) = Jde^^ its Hilbert transform. 

Introducing the inverse function such that D[R{g)] = g, we can invert the relation above 
to obtain the hybridisation function as a functional of the local G for ?7 = 0: 

Aa[i(0„;G] = G~^ +Ao[G] - aR [aG] (98) 
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So that the lattice Green's function is also expressed as a functional of G as: 

Ga(KiO)n) = 

Inserting this into (87), we can evaluate the kinetic energy: 

E^^G„(k) = '-jde^^^ = i [-1 + aGR(aG)] 

and hence the DMFT approximation to ^[G] : 

fl 



(99) 



(100) 



DMFT 



[G]= I 
Jo 



daTr 



G{i(On)R{aG{i(0„)) 



(101) 



So that the total functional reads, in the DMFT approximation: 



TdmpAG] = Fi,„p[Ao[G]]-Tr (GAo[G]) + 

■If r 

G{i(On)R{aG{io)n)) 

oc 



+ / daTr 



(102) 



In the case of an infinite-connectivity Bethe lattice, corresponding to a semi-circular 
d.o.s of width 4t, one has: R[g] = t^g + 1/g, so that the result (92) is recovered from 
this general expression. I note that the DMFT approximation to the functional ^[G] is 
completely independent of the interaction strength U. 

The equilibrium condition (89) 5T/ 5G = thus reads , in the DMFT approxima- 
tion [3]: 

Ao[i(On, G] \dmft = R[G{ia)n)] - (103) 



G{i(On) 



This can be rewritten in a more familiar form, using (96) 

D{e) 



G{i(iin)= j 



\ 7 r ^ ' with: L,>„p = ^0-1 -G-i 



(104) 



The self-consistency condition is equivalent to the condition Acj^i[G] = 0, as expected 
from the fact that Aa=i — dY/dG. Hence, within the DMFT approximation, the lattice 
Green's function is obtained by setting a = 1 into (95): 



G(k, i(On) \dmft = 



(105) 



3.4. The Baym-Kadanoff viewpoint 

Finally, let me briefly mention that the DMFT approximation can also be formulated 
using the more familiar Baym-Kadanoff functional. In contrast to the previous section. 



When deriving this equation, it is useful to note that R{aG) + aGR'{aG) 
Dynamical Mean-Field Theory February 2, 2008 



= da[aR{aG)]. 



35 



this is a functional of all components Gij of the lattice Green's function, not only of the 
local one Gu. The Baym-Kadanoff functional is defined as: 

aBK[Gij,'Lij] = -tr In [{i(On + ^)dij - tij - Lz/io^)] - tr [L • G] + ^Lw[{Gij]] (106) 

Variation with respect to E,j yields the usual Dyson's equation relating the Green's func- 
tion and the self-energy. The Luttinger-Ward functional ^jy/ has a simple diagrammatic 
definition as the sum of all skeleton diagrams in the free-energy. Variation with respect 
to Gij express the self-energy as a total derivative of this functional: 

5Gij{i(0n) 



= TTwTTTT (107) 



The DMFT approximation amounts to approximate the Luttinger-Ward functional by a 
functional which is the sum of that of independent atoms, retaining only the dependence 
over the local Green's function, namely: 

= Y,^imp[Gii] (108) 
An obvious consequence is that the self-energy is site-diagonal: 

Zij{i(On) = diji:{i(On) (109) 

Eliminating Zu amounts to do a Legendre transformation with respect to Gu, and therfore 
leads to a different expression of the local DMFT functional introduced in the previous 
section [67]: 



rDMFT[Gii\ = -trln 



{l(0n + H--^)5ij-tij 



-tr[-— - Gii]+Y,^i^p[Gii\ 



(110) 

The Baym-Kadanoff formalism is useful for total energy calculations, and will be used 
in Sec. 5.4. 



4. THE MOTT METAL-INSULATOR TRANSITION 
4.1. Materials on the verge of the Mott transition 

Interactions between electrons can be responsible for the insulating character of a 
material, as realized early on by Mott [1, 2]. The Mott mechanism plays a key role in 
the physics of strongly correlated electron materials. Outstanding examples [2, 11] are 
transition-metal oxides (e.g superconducting cuprates), fuUerene compounds, as well as 
organic conductors'"^. Fig. 6 illustrates this in the case of transition metal oxides with 
perovskite structure ABO3 [74]. 



The Mott phenomenon may also be partly responsible for the localization of f-electrons in some rare 
earth and actinides metals, see [69, 70, 71, 72, 73] and [16, 21] for recent reviews. 
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FIGURE 6. This diagram (due to A.Fujimori [74], see also [11]) can be viewed as a map of the vast 
territory of transition- metal compounds with perovskite structure ABO3. Varying the transition metal ion 
B corresponds to gradual filling of the 3d-shell. Different substitutions on the A-site can be made (A= 
Sr,Ca and A=La,Y are mainly considered in this diagram). This allows to change either the valence of the 
transition metal ion (doping), or the structural parameters in an isoelectronic manner The shaded region 
corresponds to insulating compounds, while the unshaded one corresponds to metals. This illustrates the 
key role of the Mott phenomenon in the physics of transition-metal oxides. 



A limited number of materials are poised right on the verge of this electronic instabil- 
ity. This is the case, for example, of V2O3, NiS2-xSe;c and of quasi two-dimensional 
organic conductors of the fC-BEDT family. These materials are particularly interest- 
ing for the fundamental investigation of the IVIott transition, since they offer the pos- 
sibility of going from one phase to the other by varying some external parameter (e.g 
chemical composition,temperature, pressure,...). Varying external pressure is definitely 
a tool of choice since it allows to sweep continuously from the insulating phase to the 
metallic phase (and back). The phase diagrams of (Vi-^ Crx)203 and of K"-(BEDT- 
TTF)2Cu[N(CN)2]Cl under pressure are displayed in Fig. 7. There is a great similarity 
between the high-temperature part of the phase diagrams of these materials, despite very 
different energy scales. At low-pressure they are paramagnetic IVIott insulators, which 
are turned into metals as pressure is increased. Above a critical temperature Tc (of order 
~ 450K for the oxide compound and ~ 40K for the organic one), this corresponds to a 
smooth crossover. In contrast, for T < Tc a first-order transition is observed, with a dis- 
continuity of all physical observables (e.g resistivity). The first order transition line ends 
in a second order critical endpoint at (Tc^Pc). We observe that in both cases, the critical 
temperature is a very small fraction of the bare electronic energy scales (for V2O3 the 
half-bandwidth is of order 0.5 — 1 eV, while it is of order 2000 K for the organics). 

There are also some common features between the low-temperature part of the phase 
diagram of these compounds, such as the fact that the paramagnetic Mott insulator orders 
into an antiferromagnet as temperature is lowered. However, there are also striking 
differences: the metallic phase has a superconducting instability for the organics, while 
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FIGURE 7. Left; Phase diagram of (V i _ ^ Crv)203 as a function of either Cr-concentration x or pressure 
(after[75]). Increasing x hy 1% produces similar effects than decreasing pressure by ^ 4kbar, for this 
material. Right: Phase diagram of k--(BEDT-TTF)2Cu[N(CN)2]C1 as a function of pressure (after [76]). 



this is not the case for V2O3. Also, the magnetic transition is only superficially similar : 
in the case of V2O3 , it is widely believed to be accompanied (or even triggered) by 
orbital ordering[77] (in contrast to NiS2-A:Se;i;[78]), and as a result the transition is first- 
order. In general, there is a higher degree of universality associated with the vicinity of 
the Mott critical endpoint than in the low-temperature region, in which long-range order 
takes place in a material- specific manner. 

Mott localization into a paramagnetic insulator implies a high spin entropy, which 
must therefore be quenched in some way as temperature is lowered. An obvious possi- 
bility is magnetic ordering, as in these two materials. In fact, a Mott transition between a 
paramagnetic Mott insulator and a metallic phase is only observed in those compounds 
where magnetism is sufficiently frustrated so that the transition is not preempted by 
magnetic ordering. This is indeed the case in both compounds discussed here: V2O3 
has competing ferromagnetic and antiferromagnetic exchange constants, while the two- 
dimensional layers in the organics have a triangular structure. Another possibility is 
that the entropy is quenched through a Peierls instability (dimerization), in which case 
the Mott insulator can remain paramagnetic (this is the case, for example, of VO2). 
Whether it is possible to stabilize a paramagnetic Mott insulator down to T = without 
breaking spin or translational symmetries is a fascinating problem, both theoretically 
and from the materials point of view (for a recent review on resonating valence bond 
phases in frustrated quantum magnets, see e.g [79] and [80]). The compound ^-(BEDT- 
TTF)2Cu2(CN)3 may offer [81] a realization of such a spin-liquid state (presumably 
through a combination of strong frustration and strong charge fluctuations [82]), but this 
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behaviour is certainly more the exception than the rule. 



4.2. Dynamical mean-field theory of the M ott transition 

Over the last decade, a detailed theory of the strongly correlated metallic state, and 
of the Mott transition itself has emerged, based on the dynamical mean-field theory 
(DMFT). We refer to [3] for a review and an extensive list of original references [26, 29, 
30, 83, 84, 85] We now review some key features of this theory. 

Quasiparticle coherence scale. In the metallic state, Fermi-liquid theory applies be- 
low a low energy scale e^, which can be interpreted as the coherence-scale for quasipar- 
ticles (i.e long-lived quasiparticles exist only for energies and temperature smaller than 
ep). This low-energy coherence scale is given by ep ~ ZD (with D the half-bandwith, 
also equal to the Fermi energy of the non-interacting system at half-filling) where Z is 
the quasiparticle weight. In the strongly correlated metal close to the transition, Z <^ 1, 
so that Ep is strongly reduced as compared to the bare Fermi energy. 

Three peaks in the d.o.s: Hubbard bands and quasiparticles. In addition to low- 
energy quasiparticles (carrying a fraction Z of the spectral weight), the one-particle 
spectrum of the strongly correlated metal contains high-energy excitations carrying a 
spectral weight 1 — Z. These are associated to the atomic-like transitions corresponding 
to the addition or removal of one electron on an atomic site, which broaden into Hubbard 
bands in the solid. As a result, the k-integrated spectral function A{(o) — £ijA(k, co) 
(density of states d.o.s) of the strongly correlated metal is predicted [26] to display a 
three-peak structure, made of a quasiparticle band close to the Fermi energy surrounded 
by lower and upper Hubbard bands (Fig. 8 and inset of Fig. 14). The quasiparticle part 
of the d.o.s has a reduced width of order ZD ~ ep. The lower and upper Hubbard bands 
are separated by an energy scale A. 



The insulating phase: local moments, magnetism and frustration. At strong enough 
coupling (see below), the paramagnetic solution of the DMFT equations is a Mott 
insulator, with a gap A in the one-particle spectrum. This phase is characterized by 
unscreened local moments, associated with a Curie law for the local susceptibility 
T^qXq °^ ^/T, and an extensive entropy. Note however that the uniform susceptiblity 
Xq=o is finite, of order l/Jr^U/D^. As temperature is lowered, these local moments 
order into an antiferromagnetic phase [27, 83]. The Neel temperature is however strongly 
dependent on frustration [3] (e.g on the ratio t' /t between the next nearest-neighbour 
and nearest-neighbour hoppings) and can be made vanishingly small for fully frustrated 
models. 

Separation of energy scales, spinodals and transition line. Within DMFT, a separa- 
tion of energy scales holds close to the Mott transition. The mean-field solution corre- 
sponding to the paramagnetic metal at T = disappears at a critical coupling Uc2- At this 
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FIGURE 8. Local spectral function for several values of the interaction strength in DMFT. These results 
have been obtained using the IPT approximation, for the half-filled Hubbard model with a semi-circular 
d.o.s (from Ref. [3]). Close to the transition, the separation of scales between the quasiparticle coherence 
energy (e|.) and the distance between Hubbard bands (A) is clearly seen. 



point, the quasiparticle weight vanishes (Z oc \ —U jUci) as in Brinkman-Rice theory 
On the other hand, a mean-field insulating solution is found for U > Ud, with the Mott 
gap A opening up at this critical coupling (Mott-Hubbard transition). As a result, A is a 
finite energy scale for U = Ud and the quasiparticle peak in the d.o.s is well separated 
from the Hubbard bands in the strongly correlated metal. 

These two critical couplings extend at finite temperature into two spinodal lines 
Uci{T) and UciiT), which delimit a region of the {U /D,T /D) parameter space in 
which two mean-field solutions (insulating and metallic) are found (Fig. 9). Hence, 
within DMFT, a first-order Mott transition occurs at finite temperature even in a purely 
electronic model. The corresponding critical temperature Tf is of order ~ Afi/AS, 
with AE and A5 ~ \n{lS+ 1) the energy and entropy differences between the metal 
and the insulator. Because the energy difference is small (A£' ~ [Ud — /D), the 
critical temperature is much lower than D and Uc (by almost two orders of magnitude). 



' Since the self-energy only depends on frequency within DMFT, this also implies that quasiparticles 
become heavy close to the transition, with m*/m = 1/Z. In real materials, we expect however that 
magnetic exchange will quench out the spin entropy associated with local moments, resulting in a 
saturation of the effective mass close to the Mott transition. In the regime where Ep <C J, the effective 
mass is then expected to be of order J, as found e.g in slave-boson theories. Describing this effect requires 
extensions of DMFT in order to deal with short-range spatial correlations 

Dynamical Mean-Field Theory February 2, 2008 40 




FIGURE 9. Paramagnetic phases of the Hubbard model within DMFT, displaying schematically the 
spinodal lines of the Mott insulating and metallic mean-field solutions (dashed), the first-order transition 
line (plain) and the critical endpoint. The shaded crossover lines separating the different transport regimes 
discussed in Sec. 3 are also shown. The Fermi-liquid to "bad metal" crossover line corresponds to the 
quasiparticle coherence scale and is a continuation of the spinodal UdiT) above Tc. The crossover into 
the insulating state corresponds to the continuation of the Ud spinodal. Magnetic phases are not displayed 
and depend on the degree of frustration. Figure from Refs. [86] and [87]. 



Indeed, in V2O3 as well as in the organics, the critical temperature corresponding to the 
endpoint of the first-order Mott transition line is a factor of 50 to 100 smaller than the 
bare electronic bandwith. 



4.3. Physical properties of the correlated metallic state: DMFT 

confronts experiments 

4.3.1. Three peaks: evidence from photoemission 

In Fig. 10, we reproduce the early photoemission spectra of some transition metal 
oxides, from the pioneering work of Fujimori and coworkers [88]. This work established 
experimentally, more than ten years ago, the existence of well-formed (lower) Hubbard 
bands in correlated metals, in addition to low-energy quasiparticles. This experimental 
study and the theoretical prediction of a 3-peak structure from DIVIFT [26] came in- 
dependently around the same time. However, back in 1992, the existence of a narrow 
quasiparticle peak in A{(o) resembling the DIVIFT results was, to say the least, not ob- 
vious from these early data. Further studies [89] on Cai_xSrTV03 therefore aimed at 
studying the dependence of low-energy quasiparticle spectral features upon the degree 
of correlations. One of the main difficulty raised by these photoemission results is that 
the weight Z of the low-energy quasiparticle peak estimated from these early data is 
quite small (particularly for CaVOs), while specific heat measurements do not reveal 
a dramatic mass enhancement. This triggered some discussion [89, 90, 11] about the 
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Binding Energy (eV) 

FIGURE 10. Photoemission spectra of several transition metal oxides, reproduced from Ref. [88]. 
The effects of correlations increases from ReOa (a weakly correlated metal) to YTiOa (a Mott insulator). 
The plain lines are the d.o.s obtained from band structure calculations. A lower Hubbard band around 
— 1.5 eV is clearly visible in the most correlated materials, both in the metallic and insulating case. 



possibility of a strong ^-dependence of the self-energy. A decisive insight into this ques- 
tion came from further experimental developments by Maiti and coworkers [91, 18] in 
which it was demonstrated that the photoemission spectra are actually quite sensitive to 
the photon energy. Studies at different photon energies allowed these authors to extract 
the estimated spectra corresponding to the bulk and the surface of the material. Surface 
and bulk spectra were found to be very different indeed: the surface of CaVOa being 
apparently insulating-like while the bulk spectrum did show a much more pronounced 
quasiparticle peak. Very recently, high resolution, high-photon energy photoemission 
studies [92, 93] clarified considerably this issue. The high photon-energy spectrum re- 
produced on Fig. 1 1 displays a clear quasiparticle d.o.s at low-energy (with a weight 
in good agreement with m/m* and a height comparable to the LDA d.o.s), as well as a 
lower Hubbard band carrying the rest of the spectral weight. Moreover, recent calcula- 
tions [94, 95, 93] combining electronic structure methods and DMFT (see next section) 
compare favorably to the experimental spectra, on a quantitative level. 
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FIGURE 11. (a) Bulk V 3d spectral functions of SrVOs (closed circles), Sro.5Cao.5VO3 (solid line) and 
CaV03 (open squares), (b) Comparison of the experimentally obtained bulk V 3d spectral function of 
SrVOs (closed circles) to the V 3d partial density of states for SrVOs (dashed curve) obtained from the 
band-structure calculation, which has been broadened by the experimental resolution of 140 meV. The 
soUd curve shows the same V 3d partial density of states but the energy is scaled down by a factor of 0.6. 
Figure and caption from Ref. [92] (see also [93]). 

In the case of NiS2-xSe;c, angular resolved photoemission have revealed a clear 
quasiparticle peak, with strong spectral weight redistributions as a function of tempera- 
ture [96]. For the metallic phase of V2O3 , high photon energy photoemission proved to 
be an essential tool in the recent experimental finding of the quasiparticle peak (Fig. 12) 
by Mo et al. [97]. 



4.3.2. Spectral weight transfers 

The quasiparticle peak in the d.o.s is characterized by an extreme sensitivity to 
changes of temperature, as shown in the inset of Fig. 14. Its height is strongly reduced as 
T is increased, and the peak disappears altogether as T reaches e^, leaving a pseudogap 
at the Fermi energy. Indeed, above e^, long-lived coherent quasiparticles no longer 
exist. The corresponding spectral weight is redistributed over a very large range of 
energies, of order U (hence much larger than temperature itself). This is reminiscent 
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FIGURE 12. Photoemission spectra of V2O3 , for various photon energies, from Ref. [97]. The highest 
photon energy spectrum, corresponding to the greatest bulk sensitivity, reveals a prominent quasiparticle 
peak. 



of Kondo systems [98], and indeed DMFT establishes a fomial and physical connection 
[26] between a metal close to the Mott transition and the Kondo problem. The local 
moment present at short time-scales is screened through a self -consistent Kondo process 
involving the low-energy part of the (single- component) electronic fluid itself. 

These spectral weight transfers and redistributions are a distinctive feature of strongly 
correlated systems. As already mentioned, they have been observed in the photoemission 
spectra of NiS2-xSex . They are also commonly observed in optical spectroscopy of 
correlated materials, as shown on Fig. 13 for metallic V2O3 [99] and the fC-BEDT 
organics [100]. DMFT calculations give a good description of the optical spectral weight 
transfers for these materials, at least on a qualitative level [99, 101]. 



4.3.3. Transport regimes and crossovers 

The disappearance of coherent quasiparticles, and associated spectral weight trans- 
fers, results in three distinct transport regimes [99, 102, 101, 103, 87] for a correlated 
metal close to the Mott transition, within DMFT (Figs. 9 and 14): 

• In the Fermi-liquid regime T <^ £p, the resistivity obeys a law with an enhanced 
prefactor: p = pM (T/ep)-^. In this expression, pM is the Mott-Ioffe-Regel resistiv- 
ity pM °^ ha/e^ corresponding to a mean-free path of the order of a single lattice 
spacing in a Drude picture. 

• For T ~ ep, an "incoherent" (or "bad") metal regime is entered. The quasiparticle 
lifetime shortens dramatically, and the quasiparticle peak is strongly suppressed 
(but still present). In this regime, the resistivity is metallic-like (i.e increases with 
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FIGURE 13. Left: Optical conductivity of metalUc V2O3 [99] at T = llOK (thick line) and T = 

300K (thin line)). The inset contains the difference of the two spectra Aa{(o) = (7i7o^:(g)) — Go,ook{co). 
Diamonds indicate the measured dc conductivity Gdc- Dotted lines are for the insulating compounds 
V2-y03 with y = .013 at lOK (upper) and >■ = at 70K (lower). Right: Optical conductivity of K--(BEDT- 
TTF)2Cu[N(CN)2]Br at ambiant pressure [100], for T = 25K and T = 50K. For both materials, transfer 
of spectral weight from high energies to the Drude peak is clearly visible as temperature is lowered. 



T) but reaches values considerably larger than the ]VIott "limit" Pm- A Drude 
description is no longer applicable in this regime. 

• Finally, for ep <^T <^ A, quasiparticles are gone altogether and the d.o.s displays 
a pseudogap associated with the scale A and filled with thermal excitations. This 
yields an insulating-like regime of transport, with the resistivity decreasing upon 
heating (dp/dT < 0). At very low temperature, the resistivity follows an activated 
behaviour, but deviations from a pure activation law are observed at higher tem- 
perature (these two regimes are depicted as the "insulating" and "semi-conducting" 
ones on Fig. 9). 

These three regimes, and the overall temperature dependence of the resistivity obtained 
within DMFT are illustrated by Fig. 14. A distinctive feature is the resistivity maximum, 
which occurs close to the Mott transition. This behaviour is indeed observed experi- 
mentally in both Cr-doped V2O3 and the organics. In the latter case, the transport data 
obtained recently in the Orsay group are depicted on Fig. 15, and compared to DMFT 
model calculations [103, 87]. 



Within DMFT, the conductivity can be simply obtained from a calculation of 

the one-particle self-energy since vertex corrections are absent [104, 3]. However, 
a precise determination of both the real and imaginary part of the real-frequency 
self-energy is required. This is a challenge for most "impurity solvers". In prac- 
tice, early calculations [102, 99, 101] used the iterated perturbation theory (IPT) 
approximation[26]. The results displayed in Fig. 14 have been obtained with this 
technique, and the overall shape of the resistivity curves are qualitatively reasonable. 
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FIGURE 14. Left: Resistivity in the metallic phase close to the Mott transition {U ~ 2 AD), as 
a function of temperature, calculated within DMFT using the IPT approximation. For three selected 
temperatures, corresponding to the three regimes discussed in the text, the corresponding spectral density 
is displayed in the inset. Right: IPT results for the resistivity for values of U in the metallic regime 
(lower curves), the coexistence region (bold curve) and the insulating regime (upper two curves). From 
Ref. [86, 87]. 




T (K) P (bar) 



FIGURE 15. Left: Temperature-dependence of the resistivity at different pressures, for ({"-(BEDT- 
TTF)2Cu[N(CN)2]Cl. The data (circles) are compared to a DMFT-NRG calculation (diamonds), with 
a pressure dependence of the bandwidth as indicated. The measured residual resistivity po has been added 
to the theoretical curves. Right: Transport regimes and crossovers for this compound. Figures reproduced 
from Limelette ef a/. [103]. 



However, the IPT approximation does a poor job on the quasiparticle lifetime in the 
low-temperature regime, as shown on Fig. 16. Indeed, we expect on general grounds 
that, close to the transition, DImE becomes a scaling function [105] of (o/Ep and 
T/£*p, so that for T < e^^ it behaves as: lmL{G} = Qi)oc D{T/£*f oc T^/{Z^D) which 
leads to an enhancement of the coefficient of the resistivity by 1/Z^ as mentioned 
above. The IPT approximation does not capture this enhancement and yields the in- 
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FIGURE 16. Comparison between the IPT (dashed hnes) and NRG methods (plain lines), reproduced 
from Ref. [87]. The low-frequency behaviour of the inverse lifetime ImE clearly displays a critically 
enhanced curvature, which is not reproduced by IPT. 

correct result ImE//>r((» = 0) oc U'^T^/D^, as illustrated For this reason, the numerical 
renormalization group (NRG) has been used recently [103, 87] in order to perform 
accurate transport calculations within DMFT. This method is very appropriate in this 
context, since it is highly accurate at low energies and yields real-frequency data[32]. 
DMFT-NRG calculations compare favorably to transport data on organics, as shown on 
Fig. 15. 

The crossovers described here in electrical transport also have consequences for 
thermal transport. The thermopower, in particular, displays a saturation in the incoherent 
metal regime [ 1 06, 101]. This is presumably relevant for the cobalt- based thermoelectric 
oxides such as NaxCo02 . Finally, let me emphasize that an interesting experimental 
investigation of the correlations between transport crossovers (both ab-plane and c-axis) 
and the loss of quasiparticle coherence observed in photoemission has been performed 
by Valla et al. [107] for several layered materials. This study raises intriguing questions 
in connection with DMFT, and particularly its k-dependent extensions. 



4.4. Critical behaviour: a liquid-gas transition 

Progress has been made recently in identifying the critical behaviour at the Mott 
critical endpoint, both from a theoretical and experimental standpoint. It was been 
pointed early on by Castellani et a/.[108] (see also [109]) that an analogy exists with the 
liquid-gas transition in a classical fluid. This is based on a qualitative picture illustrated 
on Fig. 17. The Mott insulating phase has few double occupancies (or holes) and 
corresponds to a low-density "gas", while the metallic phase corresponds to a high- 
density "liquid" with many double-occupancies and holes (so that the electrons can be 
itinerant). 




Dynamical Mean-Field Theory 



February 2, 2008 



47 



I I I I I I I 



— ^ — ^ — ^ — ^ — 

ill I I II 

T — ^— ^ — ^ — 



4-i-*-6-H'-4-i 

I I I < I I I 

I ■ I I I II 



I I 

I I I I I I I 



FIGURE 17. Cartoon of a typical real-space configuration of electrons in the Mott insulator (left) and 
metallic (right) phase. The insulator has few double-occupancies or holes, and corresponds to a gas of 
these excitations. Fluctuating local moments exist in this phase. The metal has many double-occupancies 
and holes, corresponding to a dense "liquid". Electrons are itinerant in the metallic phase, and the local 
moments are quenched. Within DMFT this quenching is akin to a (self-consistent) local Kondo effect. 



Recently, this analogy has been given firm theoretical foundations within the framework 
of a Landau theory [68, 110, 111] derived from DMFT by Kotliar and coworkers. In 
this framework, a scalar order parameter ^ is associated with the low-energy electronic 
degrees of freedom which build up the quasiparticle resonance in the strongly correlated 
metallic phase close to the transition. This order parameter couples to the singular part of 
the double occupancy (hence providing a connection to the qualitative picture above), as 
well as to other observables such as the Drude weight or the dc-conductivity. Because of 
the scalar nature of the order parameter, the transition falls in the Ising universality class. 
In Table 1, the correspondence between the Ising model quantities, and the physical 
observables of the liquid-gas transition and of the Mott metal-insulator transition is 
summarized. 



In Fig. 18, the dc-conductivity obtained from DMFT in the half-filled Hubbard model 
(using IPT) is plotted as a function of the half-bandwith D, for several different temper- 
atures. The curves qualitatively resemble those of the Ising model order parameter as 
a function of magnetic field (in fact, D — Dc is a linear combination of the field h and 
of the mass term r in the Ising model field theory). Close to the critical point, scaling 
implies that the whole data set can be mapped onto a universal form of the equation of 
state: 



In this expression, 7 and 8 are critical exponents associated with the order parameter and 
susceptibility, respectively: (0) ~ h^/^ at T = Tc and % = d{^)/dh ~ |r — Td^'^. f± are 




(111) 
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TABLE 3. Liquid-gas description of the Mott critical endpoint. The associated 
Landau free-energy density reads r^^ + u(j)'^ — (a possible can be eliminated 
by an appropriate change of variables and a shift of (p). 



Hubbard model 


Mott MIT 


Liquid-gas 


Ising model 


D-D, 


P-Pc 


P - Pc 


Field h 
(w/ some admixture of r) 


T-T, 


T - Tc 


T - Tc 


Distance to 
critical point r 
(w/ some admixture of h) 


Low -CO 
spectral weight 


Low -co 
spectral weight 


Vg - VL 


Order parameter 
(scalar field 0) 




FIGURE 18. IPX calculation of the dc-conductivity as a function of the half-bandwith for the half- 
filled Hubbard model within DMFT, for several different temperatures. Increasing D drives the system 
more metallic. The curve at T = Tc displays a singularity (vertical slope: dot), analogous to the non- 
linear dependence of the order parameter upon the magnetic field at a second-order magnetic transition. 
Hysteretic behaviour is found for T < T, . 



universal scaling functions associated with T > Tc (resp. T < Tc). A quantitative study 
of the critical behaviour of the double occupancy within DMFT was made in Ref. [110], 
with the expected mean-field values of the exponents 7=1,5 = 3. 



Precise experimental studies of the critical behaviour at the Mott critical endpoint 
have been performed very recently, using a variable pressure technique, for Cr-doped 
V2O3 by Limelette et al. [112] (Fig. 19) and also for the K"-BEDT organic compounds 
by Kagawa et al. [113]. These studies provide the first experimental demonstration of 
the liquid-gas critical behaviour associated with the Mott critical endpoint, including a 
a full scaling [1 12] onto the universal equation of state (111). 
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FIGURE 19. Conductivity of Cr-doped V2O3 , at the critical endpoint T = T^, measured as a function 
of pressure P/Pc (Limelette et al. [112]). A characteristic sigmoidal form is found, which is well fitted by 
(7 — <7c ~ \P — Pc\^^^ (plain line). Inset: log-log scale. See Ref. [112] for a full experimental study of the 
critical behaviour, including scaling onto the universal equation of state. 




4.5. Coupling to lattice degrees of freedom 

Lattice degrees of freedom do play a role at the IVlott transition in real materials, 
e.g the lattice spacing changes discontinuously through the first-order transition line in 
(Vi_x Crx)203, as displayed in Fig. 20. In the metallic phase, the d-electrons participate 
in the cohesion of the solid, hence leading to a smaller lattice spacing than in the 
insulating phase. 



Both the electronic degrees of freedom and the ionic positions must be retained in order 
to describe these effects. In Ref. [1 14] (see also [115]), such a model was treated in the 
simplest approximation where all phonon excitations are neglected. The free energy then 
reads: 

F = \bq^^^^^+F,i[D{v)] (112) 
2 Vo 

In this expression, v is the unit-cell volume, Bq is a reference elastic modulus and the 
electronic part of the free-energy F^i depends on v through the volume-dependence of the 
bandwith. In such a model, the critical endpoint is reached when the electronic response 
function: 

Z^-^ (1.3) 

is large enough (but not infinite), and hence the critical temperature Tc of the compress- 
ible model is larger than Tf (at which x diverges in the Hubbard model). The com- 
pressibility K = (vd^F /dv^) diverges at Tc. This implies an anomalous lowering of 
the sound-velocity at the transition [116, 117], an effect that has been experimentally 
observed in the K-BEDT compounds recently [118], as shown on Fig. 21. 
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FIGURE 20. Left: Change of the lattice constant as a function of temperature for two samples of 
(Vi_;cCrx)203 with different Cr-concentrations. The discontinuous change in the lattice constant through 
the first-order transition transition line is clearly seen for x = .006, while the sample with x = .004 is 
slightly to the right of the critical point. Right: Percentage volume change of the unit-cell volume close to 
the critical line, reflecting the critical behaviour of the order parameter. Reproduced from Ref. [109] 



We emphasize that, within DMFT, an unambiguous answer is given to the "chicken 
and egg" question: is the first-order Mott transition driven by electronic or lattice degrees 
of freedom ? Within DMFT, the transition is described as an electronic one, with lattice 
degrees of freedom following up. In fact, it is aremarkable finding of DMFT that a 
purely electronic model can display a first-order Mott transition and a finite-T critical 
endpoint (associated with a diverging x)^ provided that magnetism is frustrated enough 
so that ordering does not preempt the transition. Whether this also holds for the finite- 
dimensional Hubbard model beyond DMFT is to a large extent an open question (see 
[41] for indications supporting this conclusion in the 2D case). 



4.6. The frontier: k-dependent coherence scale, cold and hot spots 

A key question, still largely open, in our theoretical understanding of the Mott transi- 
tion is the role of spatial correlations (inadequately treated by DMFT). This is essential 
in materials like cuprates, in which short-range spatial correlations play a key role (in 
particular magnetic correlations due to superexchange, leading to a strong tendency to- 
wards the formation of singlet bonds, as well as pair correlations). In the regime where 
the quasiparticle coherence scale e'p is small as compared to the (effective strength of 
the) superexchange J, the DMFT picture is certainly deeply modified. There is com- 
pelling experimental evidence that the quasiparticle coherence scale then has a strong 
variation as the momentum k is varied along the Fermi surface, leading to the formation 
of "cold spots" and "hot regions". Such effects have been found in recent studies using 
cluster extensions of the DMFT framework ([119], see also [39, 120]). 
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FIGURE 21. Relative change in the sound velocity of K--(BEDT-TTF)2Cu[N(CN)2]Cl as a function of 
temperature, at various pressures. The velocity variation is relative to the value at 90 K. Inset: position 
and amplitude of the anomaly below 230 bars. (Figure and caption from Ref. [118]). 



5. ELECTRONIC STRUCTURE AND DYNAMICAL MEAN-FIELD 

THEORY 



The possibility of using DMFT in combination with electronic structure calculation 
methods, in order to overcome some of the limitations of DFT-LDA for strongly corre- 
lated materials, was pointed out early on [3]. In the last few years, very exciting develop- 
ments have taken place, in which theorists from the electronic structure and many-body 
communities joined forces and achieved concrete implementations of DMFT within 
electronic structure calculations. The first papers [121, 122] implementing this com- 
bination appeared in 1997-1998, and the field has been extremely active since then. For 
reviews of the early developments in this field^^, see Refs. [123, 124, 125, 126]. For 
on-line material presented at recent workshops, see Refs. [127, 128, 129] 



5.1. Limitations of DFT-LDA for strongly correlated systems 

In Sec. 3.2, 1 briefly presented the basic principles of density functional theory (DFT). 
In practice, the local density approximation (LDA) to the exchange-correlation energy, 
and its extensions (such as the generalised gradient approximation) have been remark- 
ably successful at describing ground-state properties of many solids from first princi- 
ples. This is also the state of the art method for band structure calculations, with the 
additional assumption that Kohn-Sham eigenvalues can be interpreted as single-particle 
excitations. For strongly correlated materials however, DFT-LDA has severe limitations. 



This section is merely a brief introduction to the field and certainly not as an exhaustive review. 
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which we now briefly review. 

Issues about ground-state properties. Ground-state properties, such as equilibrium 
unit-cell volume, are not accurately predicted from LDA (or even GGA) for the most 
strongly correlated materials. This is particularly true of materials in which some elec- 
trons are very localized, such as the 4f electrons of rare-earth elements at ambiant or low 
pressure (Sec. 1.3.3). If these orbitals are treated as valence orbitals, the LDA leads to 
a much too itinerant character, and therefore overestimates the contribution of these or- 
bitals to the cohesive energy of the solids, hence leading to a too small unit-cell volume. 
If instead the f-orbitals are treated as core states, the equilibrium volume is then over- 
estimated (albeit closer to experimental value, in the case of rare earth), since binding 
is underestimated. Phenomena such as the volume-collapse transitions, associated with 
the partial delocalization of the f-electrons, (and associated structural changes) under 
pressure [15] are simply out of reach of standard methods. At high pressures however, 
the f-electrons recover itinerant character and DFT-LDA(GGA) does better, as expected. 
In some particular cases, the electrons are just on the verge of the itinerant/localized be- 
haviour. In such cases, standard electronic structure methods perform very poorly. A 
spectacular example is the 5 -phase of metallic plutonium in which the unit-cell volume 
is underestimated (compared to the experimental value) by as much as 35% by stan- 
dard electronic structure methods (Fig. 24) ! All these examples illustrate the need of a 
method which is able to handle intermediate situations between fully localized and fully 
itinerant electrons. I emphasize that this issue may depend crucially on energy scales, 
with localized character most pronounced at high-energy (short time) scales, and itiner- 
ant quasiparticles forming at low-energy (long time scales). 

Excitation spectra. Even though the Kohn-Sham eigenvalues and wavefunctions 
are, strictly speaking, auxiliary quantities in the DFT formalism used to represent the 
local density, they are commonly interpreted as energy bands in electronic structure 
calculations. This is very successful in many solids, but does fail badly in strongly 
correlated ones. The most spectacular difficulty is that Mott insulators are found to have 
metallic Kohn-Sham spectra. This is documented, e.g by Fig. 23, in which the LDA 
density of states of two Mott insulators, LaTiOs and YTiOs are shown. I emphasize 
that, in both compounds (as well as in many other Mott insulators), the Mott insulating 
gap has nothing to do with the magnetic ordering in the ground-state. Even though 
magnetic long-range order is found at low-enough temperatures in both materials (below 
Tf^ ~ 140 K in LaTiOa and Tc ~ 30 K in YTiOa ), the insulating behaviour and Mott 
gap (~ 1 eV for YTiOs ) are maintained well above the ordering temperature. In other 
cases (such as VO2 ), the insulating phase is a paramagnet and the LDA spectrum is 
again metallic. 

In strongly correlated metals, e.g close to Mott insulators, the LDA bandstructure is 

also in disagreement with experimental observations. The two main discrepancies are the 
following, (i) LDA single-particle bands are generally too broad. Correlation effects lead 
to band-narrowing, corresponding to a (Brinkman-Rice) enhancement of the effective 
masses of quasiparticles. This becomes dramatic in f-electron materials, where the large 
effective mass is due to the Kondo effect, a many-body process which is beyond the 
reach of single-particle theories, (ii) The spectral weight Z associated with quasiparticles 
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is reduced by correlations, and the corresponding missing spectral weight 1 — Z is found 
in intermediate or high-energy incoherent excitations. In correlated metals, as well as in 
Mott insulators, lower and upper Hubbard bands are observed, which are absent in the 
LDA density of states (e.g for SrVOs and CaVOs in Fig. 1 1 and Fig 23). 

Related correlation effects are observed also for pure transition metals, such as nickel, 
in which the LDA spectrum is unable to account for: the ~ — 6eV photoemission 
satellite, and for the correct values of the occupied bandwidth and exchange splitting 
between the majority and minority band in the ferromagnetic ground-state. 



5.2. Marrying DMFT and DFT-LDA 

In this section, I briefly describe the (happy) marriage of electronic structure methods 
and dynamical mean-field theory. I first give a simple practical formulation in terms 
of a realistic many-body hamiltonian, and keep for the next section the construction of 
energy functional. 

The first issue to be discussed is the choice of the basis set for the valence electrons. 
Since DMFT emphasizes local correlations, we need a localised basis set, i.e basis 
functions which are centered on the atomic positions R in the crystal lattice. Up to now, 
most implementations have used basis sets based on linear muffin-tin orbitals [130, 131] 
(LMTOs) ;i^LR(r) = XiiT" — R) (in which L = {l^m} stands for the angular momentum 
quantum number of the valence electrons). These basis sets offer the advantage to 
carry over the physical intuition of atomic orbitals from the isolated atoms to the 
solid. In the words of their creator, O.K. Andersen, LMTO- based electronic structure 
methods are "intelligible" because they are based on a minimal and flexible basis set 
of short-range orbitals [132]. There are several possible choices of basis even within 
the LMTO method. Basically, a compromise has to be made between the degree of 
localisation and the orthogonality of the basis set. The most localised basis set (the so- 
called "screened" or a-basis) is not orthogonal and will therefore involve'^ an overlap 
matrix On' = {Xl\Xl')- Since DMFT neglects non-local correlations, they may be the 
best one to choose. However, a non-orthogonal basis set may not be simple to implement, 
for technical reasons, when using some impurity solvers (e.g QMC). Orthogonal LMTOs 
basis sets are somewhat more extended. 

Another possibility is to use basis sets made of Wannier functions. This has been little 
explored yet in combination with DMFT. Wannier functions can in fact be constructed 
starting from the LMTO formalism by using the "downfolding" procedure (the so-called 
third-generation LMTO [132, 133]). Recently, DMFT has been implemented within a 
downfolded (NMTO) Wannier basis, and successfully applied to transition metal oxides 
with non-cubic structures. Other routes to Wannier functions (such as the Marzari- 
Vanderbilt construction of maximally localised Wannier functions [134]) might be worth 
pursuing. Given a basis set, the electron creation operator at a point r in the solid can be 



^ ' In the following, we assume an orthogonal basis set to simplify the formalism. The overlap matrix can 
be easily reintroduced where it is appropriate 
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decomposed as: 

vHr) = '£xlR{r)cl^ (114) 

The decomposition of the full Green's function in the solid: G(r,r',T — t') = 
— {T\j/{r, T)\i/^(r', t')) (as well as of any other one-particle quantity) thus reads: 

G(r,r',/(0) = Y,Y.XLR{r)GLLi^-R',i(0)XL'Rir'y (115) 

RR' LL' 

The simplest combination of DMFT and electronic structure methods uses a starting 
point which is similar to that of the LDA+U approach [135, 136]. Namely, one first 
separates the valence electrons into two groups: those for which standard electronic 
structure methods are sufficient on one hand (e.g / = 5,p in an oxide or I = s, p,d in 
rare-earth compounds), and on the other hand the subset of orbitals which will feel 
strong correlations (e.g I = d or I = f). This separation refers, of course, to the specific 
choice of basis set which has been made. In the following, 1 denote the orbitals with / 
in the correlated subset by the index a = {m, o} (and b, ■ ■ •)• Let us then consider the 
one-particle hamiltonian: 

Hks = lef |A>(A| = Y.hfl>0^)clLCi^' (116) 

A kL 

obtained from solving the Kohn-Sham equations for the material under consideration. 
The Kohn-Sham potential we have in mind is, in the simplest implementation, the one 
obtained within a standard DFT-LDA (or GGA) electronic structure calculation of the 
local density. In a more sophisticated implementation, one may also correct the local 
density by correlation effects and use the associated Kohn-Sham potential (i.e modify 
the self-consistency cycle over the local density in comparison to standard LDA, see 
below). A many-body hamiltonian is then constructed as follows: 

H = Hks-Hdc + Hu (117) 

In this expression, Hu are many-body terms acting in the subset of correlated orbitals 
only. They correspond to matrix elements of the Coulomb interaction, and will in general 
involve arbitrary 2-particle terms UabcdcliclcdCc. In practice however, one often makes a 
further simplification and keep only density-density interactions (for technical reasons, 
this is always done when using QMC as a solver). To simplify notations, we shall limit 
ourselves here to this case, and use: 

^U=\Y^Y^ UabflRanRb (1 18) 

^ R aba 



^mm' — ^rnm' ) ^ILm' — ^LLn,' — ^mm' J mm! (119) 



with: 

In this expression, 7„„„/ is the Hund's coupling. For a more detailed discussion of the 
choice of the matrix of interaction parameters, see e.g Ref. [136]. 
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The "double-counting" term Hoc needs to be introduced, since the contribution of in- 
teractions between the correlated orbitals to the total energy is already partially included 
in the exchange-correlation potential. Unfortunately, it is not possible to derive this term 
explicitly, since the energy within DFT is a functional of the total electron density, which 
combines all orbitals in a non-linear manner. In practice, the most commonly used form 
of the double-counting term is (for other choices, see e.g [137]): 



^aba = ^ab 



Hoc = I V^cl^Cba 
Raab 

C/(iV-i)-/(iV--i) 



(120) 



The many-body hamiltonian (117) is then soved using the DMFT approximation. This 
means that a local self-energy matrix is assumed, which acts in the subset of correlated 
orbitals only: 

= 5r,R. ( zJU ) 
In the DMFT framework, the local Green's function in the correlated subset: 

Gah{T-T') = -{Tcl{T)cb{T')) (122) 
is represented as the Green's function of the multi-orbital impurity model: 

S=- f dx f dx'Y,cl{T)[%\b{x-x')ct{x') + l-Y.^ab f dxna{x)nb{x) 
JO JO 2^ Jo 

(123) 

The Weiss function (or alternatively the dynamical mean-field, or effective hybridisa- 
tion function Aab = + lJ.)5ab — [^Q^]ab) is determined, as before, from the self- 
consistency condition requesting that the on-site Green's function in the solid coincides 
with the impurity model Green's function. The components of the Green's function of 
the solid in the chosen basis set read: 

[G-%v{y^M = {i(0n + li)5w-hfl, + V^S-^LL'{i(0n) (124) 

In this expression, the self -energy matrix ILiu is constructed by using the components 
of the impurity self-energy L^^ = [^Q^]ab — [Gj^plab into (121). The self-consistency 
condition relating implicitly % and Gimp finally reads: 



k 



(125) 

J ab 



Note that this involves a matrix inversion at each k-point, as well as a k- summation 
over the Brillouin zone (which does not, in general, reduces to an integration over the 
band density of states, in contrast to the single-band case). Also let us emphasize that, 
even though the self-energy matrix has only components in the subspace of correlated 
orbitals, the components of the Green's function corresponding to all valence orbitals 
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' DMFT LOOP" 




FIGURE 22. DMFT combined with electronic structure calculations. Starting from a local electronic 
density p (r), the associated Kohn-Sham potential is calculated and the Kohn-Sham equations are solved. 
The Kohn-Sham hamiltonian//^(k) is expressed in a localised basis set (e.g LMTOs). A double-counting 
term is substracted to obtain the one-electron hamiltonian Hq = H^^ — H^'^ . The local self-energy matrix 
for the subset of correlated orbitals is obtained through the iteration of the DMFT loop: a multi-orbital 
impurity model for the correlated subset is solved (red arrow), containing as an input the dynamical mean- 
field (or Weiss field %). The self-energy E„/, is combined with Hq into the self-consistency condition 
Eq. (125) in order to update the Weiss field (blue arrow). At the end of the DMFT loop, the components 
of the full, k-dependent. Green's function in the local basis set can be calculated and thus also an updated 
local density p(r). This is used (dashed arrow) as a new starting density for the Kohn-Sham calculation 
until a converged local density is also reached . Alternatively, in a simplified implementation of this 
full scheme, the DFT-LDA calculation can be converged first and the corresponding Hq injected into the 
DMFT loop without attempting to update p (r). 



{s^p,d, ■ ■ ■) are modified due to the matrix inversion. Correlation effects encoded in the 
self-energy affect the local electronic density, which can be calculated from the full 
Green's function as: 

p{r)=Y, ZLk(r) G^z.'(k, t = 0") ^^(0 (126) 

k 

In a complete implementation, self-consistency over the local density should also be 
reached [73, 138]. The general structure of the combination of DMFT with electronic 
structure calculations, as well as the iterative procedure used in practice to solve the 
DMFT equations, is summarised on Fig. 22. 
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5.3. An application to oxides 

On Fig. 23, 1 show the spectral functions recently obtained in Ref. [95] for SrVOa , 
CaVOs , LaTiOs and YTiOs . These oxides have the same formal valence of the d- 
shell (d^). The single electron sits in the t2g multiplet, and the (empty) eg doublet is 
well separated in energy. They have a perovskite structure with perfect cubic symme- 
try for the first one (Fig. 2) and increasing degree of structural distortion for the three 
others (corresponding mainly to the GdFeOs-like tilting of oxygen octahedra). These 
calculations were performed in a downfolded (NMTO) basis set, including the off- 
diagonal components of the self-energy matrix. The latter are important for the com- 
pounds with the largest structural distortions. For comparison, the LDA density of states 
are shown on the same plot. For an independent DMFT calculation of the Ca/SrVOs 
compounds, see Ref. [94, 93] and Ref. [121, 139] for early calculations of the doped 
system Lai_;cSr;^Ti03 . The spectra in Fig. 23 have features which should be familiar to 
the reader at this point, namely: 

• SrVOa and CaVOs are correlated metals with lower (~ —1.5 eV) and upper 
(~ 2.5 eV) Hubbard bands, as well as a relatively moderate narrowing of the 
quasiparticle bandwith. The calculated spectra compare favorably to the recent 
photoemission experiments of Fig. 11 (see [93] for a comparison). 

• LaTiOs and YTiOa are Mott insulators, with quite different values of the Mott gap 
(~ 0.3 eV and ~ 1 eV, respectively) as observed experimentally. It was emphasized 
in [95] that the main reason for this difference is that the orbital degeneracy of the 
?2g multiplet is lifted to a greater degree in YTiOs than in LaTiOs due to the larger 
structural distortion. Indeed, reducing orbital degeneracy is known to increase the 
effect of correlations (for comparable interaction strength) [140, 141, 142, 143]. 
It was also found in Ref. [95] that both compounds develop a very pronounced 
orbital polarization, of a quite different nature in each compound (see [144] for a 
discussion of orbital ordering in these materials and [145] for a recent experimental 
investigation). 

This example, as well as several other recent studies, demonstrate that the embedding 
of DMFT within electronic structure calculations yields a powerful quantitative tool 
for understanding the rich interplay between correlation effects and material-specific 
aspects. 



5.4. Functlonals and total- energy calculations 

In order to discuss total energy calculations in the LDA+DMFT framework'*, it 
is best to use a formulation of this scheme in terms of a (free-) energy functional. 
Kotliar and Savrasov [138, 147] have introduced for this purpose a ("spectral-density- 



^° I acknowledge a collaboration with B. Amadon and S. Biermann [146] on the topic of this section. 
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FIGURE 23. LDA+DMFT spectral densities of the transition-metal oxides discussed in the text, from 
Ref. [95]. The (QMC) calculations were made at J = 770K. For comparison, the LDA d.o.s are also 
displayed (thin lines). 



") functional of both the total local electron density p(r) and the on-site Green's 
function in the correlated subset: (denoted Gat for simplicity in the follow- 

ing). Let us emphasize that these quantities are independent, since Gab is restricted 
to local components and to a subset of orbitals so that p(r) cannot be reconstructed 
from it. The functional is constructed by introducing (see Section. 3) source terms 
■^(r) = vks{^) — Vc{r) and ALab{iC0n) coupling to the operators v/'''(r)i//(r) and to 

LRXa(.^-^)¥ir,T:)¥H^'^T^')Xb{r' = CaR(T)cJ,j(T')> respectively. Furthermore, 
the Luttinger-Ward part of the functional is approximated by that of the on-site local 
many -body hamiltonian Hu — Hue introduced above. This yields: 

^[P{^)^Gab;VKs{^),^ab]LDA+DMFT = 

-tr ln[i(Dn + ju + 5 V2 - VKs{r) - X*-^ -X] -Jdr {vks - Vc)p (r) - tr [G.AL] + 
+i Jdrdr'p (r)C/(r - r')p (r') + E,,[p (r)] + Lr (^>,-.p[G««] - ^>z)c[G««]) 

In this expression, x*-^^ X denotes the "upfolding" of the local quantity AZ to the 
whole solid: X*-^^X = LRLahXai^ ~ R)Lab{iC0n)Xb{r' - R). Variations of this func- 
tional with respect to the sources 5Q./5vks — and 5D./ SLab — yield the standard 
expression of the local density and local Green's function in terms of the full Green's 
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function in the solid: 

p(r) = (r|G|r) , Gabim) = {XaR\G\XbR) (127) 



with: 

G 



1 

i(On + ii + -V^- VKs{r) - X*-^ X 



(128) 



or, in the local basis set (see (124)): 

1 

LL' 



G= Y,\XL^.)\{i(On + ^^).l-h^^iK-^'^ii(On)] ,{XL'^.\ (129) 



k,LL' 

From these relations, the Legendre multiplier functions v^-^ and AL could be eliminated 
in terms of p and Gab, so that a functional of the local observables only is obtained: 

^LDA+DMFT [p , Gab] = ^LDA+DMFT [p (r) , G^^ ; A [p , G] , AL [p , G] ] (130) 

Extremalisation of this functional with respect to p (dF/ dp — 0) and Gab (^r/ ^ Gab — 
0) yields the expression of the Kohn-Sham potential and self-energy correction at self- 
consistency: 

v;,5(r) = v,(r) + J Jr'C/(r-r')p(r') + ^ (131) 

^^'-Tg-;-Jg-^=^-^ ^^^^^ 

Hence, one recovers from this functional the defining equations of the LDA+DMFT 
combined scheme, including self -consistency over the local density (127). Using (66) 
and (61), one notes that the free-energy can be written as: 

^LDA+DMFT = ^oFr + tr In Gj^5(k, /o^)"^ - tr In G(k, /o^)-! -tr[Gi„pi:''"P] +LR^>;mp + 

+tT[Gin,pV'^]-l,^^DC (133) 

In this expression, CIdft is the usual density-functional theory expression (66), while 
Gks is the Green's function corresponding to the Kohn-Sham hamiltonian, i.e without 
the self -energy correction: 

G]^l = icOn + ^-hKsO^) (134) 

A careful examination of the zero-temperature limit of (133) leads to the following 
expression of the total energy [146]: 

Elda+dmft = Edft — El + (Hks) + (Hu) — Edc (135) 

= Edft + Lk,zx' ^fZ' [(4k<^i'k)£'MFr - {cli,CL'k)Ks] + {Hu) - Edc (136) 

The first term, Edft is the energy found within DFT(LDA), using of course the local 
density obtained at the end of the LDA+DMFT convergence cycle, namely: 

^DFr =£ef + J dr[vc{r)-VKs{r)]p{r) + ^Jdrdr'pir)u{r-r')p{r')+E,c[p] 
^ (137) 
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Hence the total energy within LDA+DMFT is made of several terms. Importantly, it 
does not simply reduce to the expectation value {H) of the many-body hamiltonian 
(1 17) introduced in the previous section. Furthermore, {Hks) — tr [HksG] must be eval- 
uated with the full Green's function including the self-energy correction. Therefore, 
this quantity does not coincide with the sum of the (occupied) Kohn-Sham eigenvalues 
Y!x = IvHksGks- Eq. (135) expresses that the latter has to be removed from Edft, in 
order to correctly take into account the change of energy coming from the Kohn-Sham 
orbitals. This change can also be written {Hxs)dmft — {Hks)ks = ^^[{G — Gks)Hks]- 
This is used in the second expression for the energy, which emphasizes the modi- 
fication of the density matrix (c|j^c/^'k) by correlations. Finally, the double-counting 
correction to the energy is the zero-temperature limit of — (Hoc) + tr [GE^''] — 4>dc- 
The simplest form of double-counting correction (neglecting / for simplicity) cor- 
responds to: ^Dc[Gab] = UN{N-\)/l with N = Zana = L«trG«a. Hence = 
d^Dc/dGat = U{N-\/l)na5ab, and {Hoc) = \x[GYP^] = UN{N-\/l) so that, fi- 
naily. Edc = UN{N -I) /2. 

Another formula for the total energy within LDA-i-DMFT has been used by Held et 
al. in their investigation of the volume collapse transition of Cerium [148, 149]. 

Total energy calculations within LDA+DMFT, with full self-consistency on the local 
density have been performed by Savrasov, Kotliar and Abrahams [73, 138, 147] for 
metallic plutonium with fee structure, corresponding to the 5-phase. The results are 
reproduced in Fig. 24, in which the total energy is plotted as a function of the unit-cell 
volume (normalised by the experimental value), for different values of the parameter 
U . It is seen that the GGA calculation underestimates the volume by more than 30%. 
As U increases, the minimum is pushed to higher volumes, and good agreement with 
experiments is reached for U in the range 3.8 — 4 eV. Interestingly, in the presence 
of correlations, the energy curve develops a metastable shallow minimum at a lower 
volume, which can be interpreted as a manifestation of the a-phase (which has a more 
complicated crystal structure however). For the corresponding spectra, see [147]. In 
these DMFT calculations, the 5-phase of plutonium is described as a paramagnetic 
metal, in agreement with experiments. In contrast, a static LDA-I-U treatment[72, 150] 
also corrects the equilibrium volume, but at the expense of introducing an unphysical 
spin polarization 



5.5. A life without U: towards ab-initio DIVIFT 

The combination of DMFT with electronic structure methods described in the previ- 
ous section introduces a matrix U of local interaction parameters acting in the subset of 
correlated orbitals, as in the LDA+U scheme. Some of these parameters can be deter- 
mined from constrained LDA calculations, or instead they can be viewed as adjustable. 



^ For an alternative description of the 5-phase of plutonium, in which a subset of the f-electrons are 
viewed as locaUsed while the others are itinerant, see [151] 
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FIGURE 24. Total energy of fee plutonium as a funetion of unit-eell volume (normalised by the 
experimental volume of the 5-phase), reproduced from Ref. [73, 147]. The upper curve is the GGA result. 
Other curves are from LDA+DMFT with different values of U . The lower curve is for the bcc structure. 



Furthermore, introducing these interactions implies the need for a "double-counting" 
correction in order to remove the contribution to the total energy already taken into ac- 
count in the (orbital-independent) exchange correlation potential. As such, this theory 
has great practical virtues. However, going beyond this framework and being able to 
treat the electron-electron interaction entirely from first- principles is a tempting and 
challenging project. Work in this direction have appeared recently [152, 125, 153, 154, 
155, 156, 6]. 

Physically, the Hubbard interaction is associated with the screened Coulomb interac- 
tion as seen by a given atom in the solid. Screening, is essential for estimating the order 
of magnitude of this parameter correctly. The naive view that IJ is simply the on-site 
matrix element of the Coulomb potential in the local basis- set would lead to values 
on the scale of tens of electron-volts, while the appropriate value in the solid is a few 
eV's ! This immediately points towards a key notion: that, in fact, the Hubbard f/ is a 
concept which depends on the energy-scale. At high energies (say, above the plasma fre- 
quency in a metal), it has a very large value associated with the bare, unscreened, matrix 
element, while at low energy screening takes place and it is considerably reduced. For 
first-principle RPA studies of the frequency dependence of the screened local interaction, 
see [5, 6]. 

In fact, the screened effective interaction in a solid can be related, quite generally, to 
the density-density correlation function. Let us start from the first-principles hamilto- 
nian: 



^ = - L- 5 V2 + 1, v(rO + \ L-^,. uixi - Xj) 
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= -ydr\i/'^'^^\i/ + Idrv{r)n{r) + ydrdr'u{r-r') : n(r)n(r') : (138) 



in which M(r — r') = e^/|r — r'| is the bare Coulomb interaction, n = v/^(r)v/(r) and 
: ( ) : denotes normal ordering. The (connected) density-density correlation function is 
defined as: 

X{r, r'; t - t') = {T (n(r, t) - p(r)) {n{r', x') - p {r')) (139) 
with p(r) = (n(r)) the local density. The screened effective interaction reads: 

W{yy ,i(o) ^u{r-r')- j dr\dY2u{Y -Yi)x{ri-r2;i(o)u{r2-r') (140) 

This can also be expressed in terms of the polarization P = —X-[^ — u.x]~^ as W = 
u.[l — P.u\~^ (the dot is an abbreviation for spatial convolutions). We emphasize that 
in this expression, P is the exact polarization operator, not its RPA approximation. 
The screened interaction W can be interpreted as the correlation function of the local 
scalar potential field conjugate to n{r), as can be shown from a Hubbard-Stratonovich 
transformation. 

Armed with this precise formal definition of the screened interaction in the solid (and, 
naturally, also of the full Green's function G(r,r';T- t') = -{T\if{r, T)\if^{r',T'))), we 
would like to adopt now a local picture in which we focus on a given atom. This is 
done,as before, by specifying a complete basis set of functions Xlr{^) localised around 
the atomic positions R. There is of course some arbitrariness in this choice, as already 
discussed. Adopting a local point of view, we focus on the matrix elements of the Green's 
function and of the screened effective interaction on a given atomic site: 

GabiiCO) = {XaR\G\XbR) , Wa,a2a3a4{i(0) = {XaiRXajRlWlXaiRXa^R) (141) 

In this expression, the indices a,b, - - can run over the full set of valence orbitals, or 
alternatively over a subset corresponding to the more strongly correlated ones. This is 
a matter of choice of the local quantities we decide to focus on. Following the point of 
view developed in the third section of these lectures, the key idea is again to introduce 
an exact representation of these local quantities as the solution of an atomic problem 
coupled to an effective bath. Because we want to represent the local components of 
both G and W , this effective problem now involves two Weiss functions, both in the 
one-particle and two-particle sectors. This is an extended form of dynamical mean-field 
theory (EDMFT). The action of the local problem reads: 

S = Jdxdx' [-Lctir)^-,'{^-^')cb{t') + 
H L : <(t)q,(t) : %,a2a,a,{T-T') : < (tOq,(tO :] (142) 

The local screened interaction is calculated from this effective action as: Wj,„p = — 
Ximp'^ with Ximp the 2-particle impurity correlation funcition. The two Weiss fields ^ 
and are adjusted in such a way that Gimp — Gab ^^d Wimp — Wabcd^ the local quantities 
in the solid. The impurity model (142) can be viewed as an atom hybridised with an 
effective bath of non-interacting fermions and also coupled to a bath of fluctuating 
electric scalar potentials. 
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This construction provides an unambiguous definition of the Hubbard interactions 
'^abcd{iOi) in the solid (as well as of the usual dynamical mean-field assuming 
of course that the local components of the screened interaction W and of the Green's 
function G are known. Frequency- dependence of ^ is essential in a proper definition 
of these Hubbard interactions, at least when a wide range of energy scale is considered. 
Naturally, one degree of arbitrariness remains, associated with the choice of the basis 
set: ^ will change when a different basis set is considered, keeping the same form of 
the effective interaction W{r,r';i(o) in the full solid. 

To proceed from these formal considerations to a practical scheme, we need to decide 
how W and G will actually be calculated, and this of course will involve approximations. 
Again, a free-energy functional is an excellent guidance and indeed such a functional 
of the full G{r,r';i(o) and W{r,r';i(o) has been introduced by Almbladh et al.[157], 
generalizing the Baym-Kadanoff construction (see also [47] for independent work). The 
functional reads: 

r{G,W) = TrlnG - Tr[{G]j^ -G-^)G] - ^TrlnW + ^Tr[{u-^ - W-'^)W] +'¥[G,W] 

(143) 

Gfj^ = icOn + /i + — vh corresponds to the Hartree Green's function with vh being 
the Hartree potential. For a derivation of (143) using a Hubbard-Stratonovich trans- 
formation and a Legendre transformation with respect to both G and W, see [47]. 
The functional ^[G,W] is a generalization of the Luttinger-Ward functional ^[G], 
whose derivative with respect to G gives the self-energy. Here we have, similarly (from 

dr/dG^dr/dw ^oy. 

XXV XUJ 

G^ = Gr,^-L''- , E^^ = ^ ; = M 1 - P , P = -2-^— (144) 
" ' dG dW 

A well established electronic structure calculation method, which offers in part an 
alternative to DFT-LDA, is the so-called GW approach [158] (see [159] for a review). 
This corresponds to the following approximation to the ^-functional: 

^GWA drdr' j dxdx' G{ry ,x-x')W{ry ,T-x')G{r' ,r,x' -x) (145) 

which yields the RPA-like approximation to the polarisation and exchange-correlation 
self-energy: P = G* G and E-^*^ = — G* W. The GW approximation to the ^-functional 
is easily written in terms of the components of G and W in the chosen basis set: 

^GWA = -\fdxl^l^ CL;(^)<E.;4(^)Gf;j,(-T) (146) 

This can be separated into a contribution ^ctka''"^ from non-local components (corre- 
sponding to the terms with R ^ R' in (146)) and a contribution ^q^^ [G^^, W^^] from 
local components only (R = R'). 

The GW approximation does treat the screened Coulomb interaction from first- 
principles, but does not treat successfully strong correlation effects. Recently, it has 
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been suggested to improve on the GWA for the local contributions by using the DMFT 
framework [152, 153] (see also [125, 154, 155, 156]). One can think of different approx- 
imations to the ^-functional in this context, depending on whether the DMFT approach 
is used for all the valence orbitals L — s,p,d,---, or for a subset (corresponding to the 
index ...) of correlated orbitals only. The corresponding ^-functional reads: 

R 

(147) 

In this expression, ^imp is the ^-functional corresponding to the local effective model 
(142), while removes the components from ^'^^^ which will be taken into account 
in ^imp, namely: 

= -^I /rfT £ G« (T)wS(T)Gg«(-T) (148) 

^ R abed 

If all valence orbitals are included in the DMFT treatment, the second term in the r.h.s of 
(147) is absent altogether. If only a correlated subset is treated with DMFT, AM^ can be 
thought of as a term preventing double-counting of interactions in the correlated subset. 
In this context however, in contrast to LDA+DMFT, the form of this double-counting 
correction is known explicitly. 

Taking derivatives of this functional with respect to the components of G and W, one 
sees that, in the GW+DMFT approach, the non-local components of the self-energy 
and of the polarization operator keep the same form as in the GWA, while the local 
components are replaced by the ones from the effective impurity model (possibly in the 
correlated subset only). The GW+DMFT theoretical framework is fully defined by (147) 
and the form of the impurity model (142). As before, an interative self -consistent process 
must be followed in order to obtain the self-energy and screened effective interaction, 
as well as the dynamical mean-field ^ and effective Hubbard interactions "^Z . This is 
described in more details in Refs. [152, 155, 156]. Concrete implementations of this 
scheme to electronic structure (and to model hamiltonians as well) is currently being 
pursued by several groups. For early results, see [152, 153, 154, 155, 156, 6]. 



6. CONCLUSION AND PERSPECTIVES 

In these lectures notes, I have tried to give an introduction to some aspects of the 
physics of strong electron correlations in solids. Naturally, only a limited number of 
topics could be covered. The field is characterized by a fascinating diversity of material- 
dependent properties. It is, to a large extent, experimentally driven, and new discoveries 
are undoubtedly yet to come. Also, new territories outside the traditional boundaries 
of solid-state physics are currently being explored, such as correlation effects in nano- 
electronic devices or the condensed matter physics of cold atoms in optical lattices. 

On the theory side, these lectures are influenced by the author's prejudice that (i) 
physics on intermediate energy scale matters and may be a key to the unusual behaviour 
of many strongly correlated materials and that (ii) quantitative theoretical techniques 
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are essential to the development of the field, in combination with phenomenological 
considerations and experimental investigations. 

Dynamical mean-field theory is a method of choice for treating these intermediate 
energy scales. The basic principles of this approach have been reviewed in these lectures. 
On the formal side, analogies with classical mean-field theory and density-functional 
theory have been emphasized, through the construction of free-energy functional of 
local observables. A distinctive aspect of DMFT is that it treats quasi-particle excitations 
and higher energy incoherent excitations, on equal footing. As a result, it is able to 
describe transfers of spectra weight between quasiparticle and incoherent features as 
temperature, coupling strength, or some other external parameter (doping, pressure,...) 
is varied. I have emphasized that the quasiparticle coherence scale plays a key role in 
the physics of a strongly correlated metal. Above this scale, which can be dramatically 
reduced by correlations, unusual (non-Drude) transport and spectroscopic properties are 
observed, corresponding to an incoherent metallic regime. This is the case, in particular, 
for metals which are close to a Mott insulating phase. I have briefly reviewed the DMFT 
description of these effects in these lectures, in comparison to experiments, as well as 
the detailed theory of the Mott transition which has been one of the early successes 
of this approach. I have also provided an (admittedly quite succinct) introduction to the 
recent combination of DMFT with electronic structure calculations. These developments 
have been made possible by researchers from two communities joining forces towards a 
common goal. It provides us with a powerful quantitative tool for investigating material- 
dependent aspects of strong electron correlations. 

Despite these successes, some key open questions in the physics of strongly correlated 
electron systems remain out of reach of the simplest version of DMFT. Indeed, in 
materials like cuprates, short-range spatial correlations play a key role (in particular 
magnetic correlations due to superexchange, leading to a strong tendency towards the 
formation of singlet bonds, as well as pair correlations). These correlations deeply 
affect the nature of quasiparticles. There is compelling experimental evidence that the 
quasiparticle coherence scale has thus a strong variation as the momentum k is varied 
along the Fermi surface, leading to the formation of "cold spots" and "hot regions". 
Extending the DMFT framework in order to take these effects into account may well be 
the most important frontier in the field. 
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